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ABSTRACT 


The  primary  objective  of  the  development  of  the  Variable  Point  Guidance  and  Targeting 
Technique  is  to  devise  a  quick  reaction  guidance  and  targeting  system  for  space  vehicles 
which  does  not  limit  the  intrinsic  propulsion  capability  of  the  booster  and  upper  stages  in 
applications  involving  arbitrary  rendezvous  missions.  Logic  and  equations  for  such  a 
system  which  operate  in  real  time  and  which  fit  a  modern  aerospace  computer  have  been 
developed  and  simulated.  The  system  is  provided  with  the  ability  of  accepting  a  revised 
target  ephemeris  or  tbj  ephemeris  of  a  completely  new  target  aftei  lift-off,  giving  the 
operational  .flexibility  which  permits  in-flight  changing  of  mission  objectives  at  any  time. 

The  purpose  of  the  present  work  has  been  (1)  to  improve  the  efficiency  and  effectiveness  of 
the  techniques,  (2)  to  verify  the  modifications  by  a  detailed  simulation  containing  a  repre¬ 
sentation  of  a  typical  booster  and  upper -stage  configuration,  (3)  document  the  complete 
program  and  (4)  to  code  an  aerospace  computer  to  the  extent  that  an  accurate  sizing  and 
timing  estimate  can  be  made.  The  program  was  exercised  by  complete  simulation  of  twelve 
rendezvous  missions  against  four  widely  different  tareet  situations.  An  evaluation  of  the 
results  of  the  rendezvous  simulation  studies  indicates  that  the  efficiency,  effectiveness  and 
general  a;  plicability  of  the  computer  program  was  improved  and  extended  so  that  any 
satellite  rendezvous  mission  thai  is  within  the  potential  capability  of  Urn  booster  and  its 
associated  hardware  can  successfully  be  initiated  and  completed  by  the  insertion  of  the  tar¬ 
get  ephemeris. 
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SECTION  T 


INTRODUCTION  AND  SUMMARY 


1.  Introduction.  Today's  capability  for  operating  in  space  requires  consideration  of  more 
ambitious  rendezvous  missions  than  those  carried  out  in  the  Gemini  program.  The  need  to 
"  able  to  rescue  astronauts  stranded  in  space  requires  the  use  of  a  self-contained  rendez¬ 
vous  guidance  technique  which  can  instantaneously  plan  and  execute  rendezvous  with  an 
arbitrary  target.  Launching  interceptor  vehicles  from  aircraft,  space  stations  or  other 
mobile  sites,  and  allowing  astronauts,  in  orbit,  -o  change  the  missile  to  encompass  pre¬ 
viously  unplanned  objectives,  also  call  for  more  flexibility  in  flight  planning  and  guidance. 

The  desirable  characteristics  o f  a  guidance  scheme  to  handle  these  and  other  generalizations 
of  the  space  rendezvous  mission  are  that  it  be  self -targeting,  have  fast  reaction  time,  be 
free  from  dependence  on  ’’nominal"  trajectories,  depend  on  a  minimum  number  of  vehicle 
characteristics,  operate  in  real  time,  accept  target  ephemeris  changes  after  launch,  and  be 
efficient,  flexible,  self-contained  and  compact. 

Self -targeting  means  having  the  capability  of  solving  the  targeting  problem  within  the  on¬ 
board  guidance  computer,  without  aid  from  a  ground-based  computer.  The  targeting  pro¬ 
blem  involves  determining  when  to  launch,  where  to  rendezvous,  and  when  and  where  to  per¬ 
form  every  intermediate  maneuver.  Fast  reaction  refers  to  the  capability  for  being  able  to 
solve  the  targeting  problem  rapidly,  say  in  one  second  or  less,  so  as  not  to  delay  the  launch. 

The  guidance  should  be  free  from  dependence  on  nominal  trajectories  so  that  it  can  handle 
any  target  whose  orbit  is  within  the  propulsion  capability  of  the  booster.  It  would  depend  on 
the  minimum  number  of  vehicle  characteristics*  so  that  the  payload  weight  or  other  vehicle 
parameters  may  be  changed  right  up  to  the  tim*.  of  launch  witnout  requiring  recomputation 
of  guidance  parameters.  It  should  operate  in  real  time  to  guide  the  vehicle  in  a  stab’e  man¬ 
ner  through  all  phases  cf  flight,  correcting  for  perturbations  that  were  not  predicted  when 
the  flight  plan  was  developed  and  modifying  the  flight  plan  if  necessary.  In  a  fast  reaction 
situation,  the  target  data  may  be  as  good  before  launch  as  after,  due  to  the  acquisition 
of  more  measurements  on  the  target  by  the  tracking  network.  Therefore,  it  is  desirable 
that  the  system  can  accept  an  updated  target  ephemeris  after  launch  and  compute  an  effective 
and  efficient  course  of  action  that  will  reliably  compute  the  mission.  The  program,  iu  be 
useful,  must  be  compact  enough  to  fit  the  available  memory  and  spec  i  capabilities  of  existing 
Aerospace  computers. 

These  capabilities  are  provided  to  a  high  degree  by  the  Variable  Point  Guidance  and  Target¬ 
ing  Program.  The  approach  has  been  to  integrate  the  functions  of  targeting,  flight  planning 
and  guidance.  The  function  of  targeting  is  to  predict  the  target  state  vector  at  any  desired 
time  and  also  find  an  accepted  flight  plan  for  possible  use.  The  function  of  flight  planning 
is  to  determine  the  best  or  most  efficient  plan  for  accomplishing  rendezvous. 
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The  flight  plans  that  are  used,  take  into  account  the  interceptor  characteristics  combined 
with  the  principal  results  of  optimization  theory  (such  as  Hchmann  transfers)  rather  than 
mechanizing  the  theory  itself.  Heuristic  reasoning  is  used  to  make  a  reasonable  tradeoff 
between  time  cf  flight  and  amount  of  propellant  consumed  over  a  wide  class  of  targets. 

The  function  of  guidance  is  to  generate  the  necessary  command  and  control  signals  so  that  the 
prescribed  flight  plan  is  followed.  The  function  of  guidance  in  the  Variable  Point  Guidance 
and  Targeting  Program  is  not  only  to  generate  the  guidance  commands  based  upon  an  efficient 
flight  plan,  but  also  to  alter  the  flight  plan  itself  when  it  is  ntr^ssary  or  more  efficient  to  do 
so.  This  is  accomplished  by  the  fact  that  the  targeting  routine,  flight  plan  generator  and 
guidance  equations  are  integrated  into  a  relatively  simple  routine  which  operates  in  reu.time 

As  new  theoretical  results  from  astrodynamics  become  available  and  as  computer  capability 
increases,  evolutionary  improvements  can  be  made  to  the  flight  plan  logic  and  the  guidance 
equations  to  approach  closer  and  closer  to  optimal  performance  in  a  still  wider  class  of 
cases. 

2.  Summary.  Previous  work  on  Variable  Point  Guidance  and  Targeting,  contract  AF04(695)- 
273  and  AJ  04(695) -633  were  involved  with  the  development  of  guidance  techniques  applicable 
to  the  general  rendezvous  mission.  The  first  contract  AF04(695)-273  was  a  follow  on  to 
SAINT  Program  706.  Variable  Point  Guidance  and  Targeting,  a  form  of  Explicit  Point  Guid¬ 
ance,  was  developed  out  of  the  recognized  need  for  more  flexibility  and  faster  reaction  time. 

It  was  observed  early  in  Program  706  in  working  with  Explicit  Point  and  other  Guidance 
methods  that  launch  windows  were  limited  not  only  by  vehicle  propulsion  capability  but  also 
by  the  limitations  inherent  in  the  targeting,  flight  planning  and  guidance  method.  This  led  to 
combining  the  targeting,  flight  planning  and  guidance  into  a  single  program  as  shown  in 
Figure  1.  In  order  to  maintain  efficiency  and  flexibility,  the  aim  points  and  rendezvous  point 
for  a  particular  lift  -off  time  were  planned  to  be  variable.  The  phasing  was  initially  accomp¬ 
lished  by  making  use  of  parking  orbits  and  bi-eiliptic  transfers.  The  angle  out  of  plane  was 
removed  at  the  first  crossing  of  the  line  of  nodes.  Relatively  good  results  were  obtained 
from  this  simple  system  especially  for  missions  involving  near  circular  orbits. 

In  contract  AF04(695)-633  means  were  studied  for  improving  the  performance  efficiency  and 
overall  effectiveness  for  the  general  rendezvous  mission.  The  plane  change  was  optimally 
partitioned  among  the  nodal  burns,  the  ability  to  change  targets  at  any  time  or  update  the 
initial  ephemeris  was  included,  the  orbit  transfer  rnd  phasing  policies  were  modified  to 
improve  operation  against  highly  elliptical  target  orbits  and  the  number  of  burn  points  was 
decreased,  therefore  increasing  reliability. 

Alternate  methods  of  rendezvous  guidance  were  studied  and  successfully  simulated  using  the 
Titan  IOC  as  a  booster.  A  report  on  these  studies  (SSD-100)  dated  15  July  1965,  Reference 
(1),  gives  the  analysis  of  the  various  methods  and  the  simulation  results  for  several  test  cases. 
The  purp  ~e  of  the  present  contract  AF04(695)-946  is  to  refine  and  document  the  Variable 
Point  Guidance  and  Targeting  equations  and  logic  as  well  as  the  simulation  program  of  the 
guided  Titan  me  multistage  booster  and  its  environment.  In  addition,  the  simulation  program 
is  to  be  implemented  with  ar>  SC4020  automatic  plotting  system  and  the  efficiency  and  effect¬ 
iveness  of  the  program  are  to  be  demonstrated  by  a  series  of  twelve  runs  against  four 
specific  targets.  The  equations  and  logic  are  also  to  be  programmed  for  sizing  in  a  modern 
aerospace  computer.  The  work  is  broken  into  the  following  four  tasks: 

a.  Task  1  -  Program  Improvement  and  Documentation.  Three  items  were  identified 
in  past  work  which  were  considered  to  be  of  primary  importance  for  improving  the  effi¬ 
ciency  in  the  general  rendezvous  mission.  The  tangent  transfer  and  full  orbit  phasing  policy 
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Figure  1.  Variable  Point  Computation  Guidance  and  Targeting 


was  studied  in  the  previous  contract  and  found  to  be  markedly  more  effective  in  missions 
involving  highly  elliptical  target  orbits.  The  subprogram  was  developed  and  required  inte¬ 
gration  into  the  system. 

The  second  item  for  improvement  involves  the  reference  vector  for  orientation  during  coast 
periods  in  preparation  for  the  next  on-orbit  burn.  The  method  previously  used  did  not 
always  produce  a  stable  reference  cr  a  correct  time  to  start  the  burn.  The  method  for 
correcting  this  defect  was  to  use  the  reference  vector  associated  with  the  burn  point  de¬ 
rived  in  the  flight  plan. 

The  third  item  concerned  the  effects  of  earth  oblateness  upon  the  targeting  and  guidance 
equations.  The  computation  for  predicting  the  future  interceptor  and  target  vehicle  position 
and  velocity  should  account  for  oblateness  effects.  This  was  done  in  the  navigation  equa¬ 
tions  and  in  the  target  prediction  equations,  however  the  correction  was  not  made  in  the 
missile  flight  plan  and  guidance  equations.  The  resulting  trajectory  therefore  would  cause 
the  missile  to  miss  the  rendezvous  point.  The  simplest  solution  to  this  problem  requires 
only  that  a  consistent  prediction  method  be  used  for  each  future  trajectory  segment.  The 
program  was  modified  to  include  Kepler  arcs  in  forward  computation  for  both  vehicles  and 
correct  navigation  and  targeting  for  past  computations.  The  relative  trajectory  errors 
consequently  tend  to  cancel  out  as  the  two  vehicles  approach  rendezvous. 

Task  1  also  includes  the  documentation  of  the  Variable  Point  Guidance  and  Targeting  Method 
and  a  programming  manual  of  the  simulation  program.  The  task  1  documentation  is 
divided  into  two  parts.  Part  1  includes  background  theory,  equations  and  logic,  derivation 
of  equations  and  simulation  results.  Part  2,  the  programmers  manual,  contains  flow 
charts,  the  Fortran  listing  of  guidance  and  targeting,  the  vehicle  and  its  environment  and 
the  implementation  of  the  SC4020  plotting  routine.  Task  1  also  includes  the  demonstration 
and  delivery  of  the  functioning  simulation  program  suitable  for  use  in  the  system  evaluation 
group  computer  facility  at  Wright- Patterson  Air  Force  Base. 

b.  Task  2  -  Implementation  of  the  SC4020  Plotting  Routine.  This  task  includes  the 
implementation  of  the  SC4020  plotting  routine  with  the  complete  guided  vehicle  simulation 
program.  In  addition,  the  complete  program  and  instructions  for  automatically  generating 
SC4020  plots  of  any  of  68  variables  are  included.  The  plot  routine  has  also  been  imple¬ 
mented  with  the  flight  plan  generator  as  an  added  feature  for  examination  of  the  launch  trade 
off  data. 


c.  Task  3  -  Generation  of  Program  Specification  for  Spaceborne  Computer  Coding. 
The  subtasks  here  include  actual  coding  of  Variable  Point  Guidance  and  Targeting  routines 
for  the  Uni  vac  1824C  aerospace  computer.  The  coding  of  the  routines  in  the  computer  is 
then  used  to  obtain  a  basis  for  estimating  the  cycle  time  and  storage  requirement  of  the 
program.  The  work  was  divided  into  the  following  subtasks: 

1.  )  Generate  the  navigation  guidance  and  targeting  equations  for 

coding  in  an  aerospace  computer 

2. )  Code  specified  Variable  Point  Guidance  and  Targeting  equations 

in  182 4C  assembly  language 

3.  )  Estimate  the  cycle  time  and  storage  requirement  of  the  Variable 

Point  Guidance  and  Targeting  Program. 

4.  )  Generate  a  document,  Volume  II  which  gives  the  results  of  the  study. 
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Using  the  results  of  the  coding,  and  making  neces.sary  comparisons  between  the  simulation 
program,  a  sizing  and  timing  study  was  also  completed.  The  purpose  was  to  verify  that 
the  Variable  Point  Guidance  and  Targeting  storage  requirements  would  fit  the  available 
space  of  an  airborne  computer.  A  timing  study  to  determine  if  all  the  calculations  could  be 
done  fast  enough  to  satisfy  the  guidance  and  targeting  requirements  was  also  completed. 

A  summary  of  the  Univac  1824  characteristics,  program  requirements  and  estimated  size 
and  time  for  the  program  are  given  in  Tables  I,  n  and  in. 

It  is  concluded  from  the  results  of  this  study  that  the  Variable  Point  Guidance  and  Target¬ 
ing  system  is  self  contained,  it  does  in  fact  fit  the  available  space  in  the  1824C  computer 
and  the  program  is  fast  enough  to  satisfy  the  guidance  and  targeting  requirements  of  a 
general  rendezvous  mission.  The  1824G  is  an  18000  word  version  of  the  1824C.  It  is  the 
same  size  but  four  pounds  heavier  with  6000  additional  memory  locations. 

d.  Task  4  -  Terminal  Guidance,  Initial  Conditions  Study.  This  task  concerns  the 
simulation  of  the  Variable  Point  Guidance  and  Targeting  system  in  specific  rendezvous 
missions.  There  are  twelve  runs  against  the  following  four  targets: 

1.  )  A  satellite  in  a  circular  orbit  at  an  altitude  of  150  nautical  miles 

inclined  at  an  angle  of  35  degrees 

2.  )  A  satellite  in  a  circular  orbit  at  an  altitude  of  1000  nautical  miles 

inclined  at  49  degrees 

3.  )  A  synchronous  satellite  in  an  elliptical  orbit  whose  perigee  is  at 

500  nautical  miles  altitude  at  an  inclination  angle  of  65  degrees 

4.  )  A  satellite  in  an  elliptical  orbit  with  altitude  of  perigee  of  200 

nautical  miles  and  altitude  of  apogee  of  300  nautical  miles  inclined 

at  an  angle  of  90  degrees. 

Three  selected  test  cases  for  each  of  the  above  targets  were  run  using  three  different  times 
of  lift-off.  The  variation  in  lift-off  time  gives  a  variation  in  the  trajectories  and  hence  a 
different  relative  velocity  at  the  beginning  of  the  terminal  phase  (when  the  terminal  sensors 
would  become  effective).  For  this  task  the  beginning  of  the  terminal  phase  was  defined  as 
the  time  that  the  relative  range  is  100  nautical  miles.  SC4020  plots  we  re  made  of  various 
parameters  of  the  boost  phase  of  one  of  the  runs,  and  plots  were  also  made  of  the  relative 
range,  relative  velocity,  line  of  sight  rate,  body  rate  commands,  specific*  force,  radius 
vector  and  velocity  vectors  for  the  last  100  nautical  miles  of  the  mission. 

Most  simulation  runs  were  well  behaved  and  terminated  within  a  reasonable  rendezvous 
criteria.  The  results  and  detailed  discussion  of  the  simulation  runs  are  given  in  Section  IV, 
Volume  I,  Part  1  of  this  report.  The  results  showed  that  modifications  in  the  terminal 
guidance  (the  present  closure  mode)  should  be  made  which  would  lower  the  sensitivity  and 
provide  alternate  means  for  carrying  out  the  flight  plan  to  assure  that  a  reliable  and  effec¬ 
tive  reference  backup  is  always  available  on  the  contingency  that  the  method  being  used 
either  fails  to  develop  an  acceptable  terminal  plan  or  indicates  that  excessive  velocity 
changes  over  those  derived  in  the  original  plan  are  required  to  complete  the  mission. 
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Table  I.  Univac  Computer  Characteristics 


Characteristics 

1824C 

1824G 

Size  -  In  Cubic  Inches 

1300 

1300 

Weight  In  Pounds 

45 

49 

Data  Word  Length  In  Bits 

24 

24 

Instruction  Length  In  Bits 

16 

16 

Memory: 

NDRO  096  (16  Bit)  Instructions 

12,096 

18,  00C 

DRO  (24  Bit)  Data  Words 

512 

512 

Storage  Cycle  Time  in  microseconds 

4 

4 

Add  Time  in  Microseconds 

8 

8 

Multiply  Time  Averages  in  Microseconds 

64 

64 

Table  II.  Storage  Requirements 


VPG  Program 

(5226  Instr  +300  Constants)  5676 

Other  (From  Flight  Plan  VII  Titan  IRC)  7858 


Total:  13,534 

Remaining  1824G:  4466 


Table  EL  Timing 


Execution  Time  (Seconds)  for  a  Typical  Burn  Cycle  0.09 
Cycle  Time  Requirement  to  be  Less  Than  2.  0 

Execution  Time  (seconds)  for  a  Typical  Coast  Cycle  C.  177 
Cycle  Time  Requirement  to  be  Less  Than  10.  0 
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3.  Conclusion  .  A  computer  program  for  the  general  solution  of  the  problem  of  guidance 
of  spacecraft  on  orbital  rendezvous  missions,  Variable  Point  Guidance  and  Targeting,  has 
been  developed,  documented,  programmed  for  sizing  on  an  aerospace  computer,  and 
demonstrated  by  simulation.  In  this  contract,  further  improvements  were  made  in  the  pro¬ 
gram  which  increased  the  efficiency  and  flexibility  and  extended  the  capability  of  the  tech¬ 
nique  to  a  wider  class  of  missions.  The  accuracy  of  the  system  has  been  improved  by  a 
better  means  to  account  for  earth  oblateness  effects.  The  more  efficient  full  orbit  phasing 
method  was  integrated  and  a  better  means  for  controlling  the  vehicle  in  relation  to  coast 
and  burn  phases  of  orbit  burns  have  been  included.  To  be  ready  for  any  early  specific  appli¬ 
cation  further  modification  and  improvements  should  be  made. 
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SECTION  n 


THEORETICAL  BACKGROUND 


1.  Introduction.  In  order  to  specify  the  optimum  trajectory  to  perform  a  given  space 
mission,  the  entire  thrust  history  (including  both  propulsive  and  aerodynamic  forces) 
must  be  determined.  This  includes  the  specification  of  the  number  of  thrusts,  their 
location  in  space  and  time,  and  the  thrust  vector  orientation  at  each  instant  during  the 
maneuver.  This  task  may  be  called  flight  plan  selection,  or  targeting,  since  it  gives  an 
overall  description  of  the  trajectory  arcs  which  comprise  the  mission.  The  velocity 
changes  needed  to  perform  the  mission  are  those  required  to  transfer  between  the 
successive  arcs  in  an  optimal  manner.  In  general,  the  velocity  requirement  is  much  more 
sensitive  to  proper  selection  of  the  flight  plan,  than  to  the  method  of  guidance  and  control 
during  the  maneuvers.  In  fact  the  fixed  axis  steering  method  used  in  the  Variable  Point 
Guidance  mode  has  been  shown  to  give  results  differing  relatively  little  from  the  optimum. 
The  central  problem  of  Variable  Point  Guidance  and  Targeting  is  therefore  the  rapid  and 
automatic  generation  of  a  near -optimum  flight  plan.  The  integration  of  efficient  methods 
of  propulsive  and  aerodynamic  maneuvers  into  the  overall  plan  constitutes  the  guidance 
part  of  the  problem. 

We  shall  here  refer  tc  targeting  as  the  means  for  providing  the  flight  plan  to  satisfy  mission 
objectives  and  to  guidance  as  that  phase  of  the  mission  which  begins  after  the  flight  plan 
has  been  chosen,  and  the  launch  process  has  begun.  During  the  guidance  phase  the  flight 
plan  is  being  revised  as  necessary,  to  take  account  of  deviations  from  the  original  plan. 

The  basic  problem  to  which  Variable  Point  Guidance  and  Targeting  has  been  addressed 
may  be  formally  stated  as  follows:  ”We  which  to  direct  a  space  vehicle  from  some  initial 
trajectory  in  spacetime  to  a  prescribed  final  trajectory,  and  to  accomplish  this  transfer 
with  minimum  AV,  minimum  time,  or  some  suitable  combination  of  the  two,  subject  to 
prescribed  constraints.”  In  this  context  ’’launch”  may  represent  an  actual  launch  from  a 
ground  site,  or  the  start  of  a  burn  from  orbit.  Similarly  "rendezvous”  may  represent  an 
actual  rendezvous  with  a  target  in  orbit,  or  achievement  of  the  desired  terminal  conditions. 
Since  AV  and  flight  time  in  general  represent  conflicting  requirements,  it  cannot  be 
established  a  priori  just  what  represents  the  most  suitable  combination  of  the  two.  In  fact, 
there  are  situations  or  missions  in  which  one  may  be  prepared  to  sacrifice  AV  to  gain 
time,  (AV  may  be  budgeted  and  time  is  to  be  minimized). 

Since  it  is  required  that  our  system  be  flexible,  it  must  anticipate  such  situations.  Thus 
the  mathematics  should  be  formulated  in  such  a  way  that  it  will  result  in  a  family  of  "ideal” 
solutions;  i.e. ,  those  for  which  rendezvous  at  various  times  is  effected  with  the  least  AV. 
These  solutions  are  then  presented  to  either  the  computer  or  a  human  operator,  who 
proceeds,  on  the  basis  of  previously  defined  logical  criteria,  to  make  the  tradeoff  be¬ 
tween  time  and  AV  which  is  appropriate  to  that  particular  mission.  Furthermore,  this 
process  should  be  performed  continually,  in  order  to  provide  a  guidance  capability  and  to 
adapt  to  situations  which  may  be  changing  in  real  time. 
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It  is  convenient  initially  to  define  a  propulsion  program  in  terms  of  V’  discrete  impulse 
burns.  There  are  then  (4n  -  6)  degrees  of  freedom.  This  may  be  seen  by  noting  that  there 
are  three  resultant  velocity  components  and  one  time  associated  with  each  impulse.  (The 
position  of  the  impulse  is  a  known  function  of  time,  as  defined  by  the  previous  trajectory.) 
This  would  give  4n  variables.  However,  six  of  these  are  constrained  when  we  prescribe 
the  space-time  target  trajectory,  three  for  target  position  at  rendezvous  time  and  three 
for  velocity;  hence  (4n  -  6)  free  s^ate  variables. 

There  are  of  course  additional  non-holonomic  constraints,  such  as  the  requirement  that 
we  not  pass  into  the  earth,  range  safety  restrictions,  etc.  However,  total  AV  may  be 
expressed  as  a  function  of  (4n  -  6)  variables. 


^L*  tR’  Pl’  P2 


P4n  -  8> 


(n-i) 


where  the  P.'s  are  independent  variables,  and  launch  time  t^  and  rendezvous  time  t„ 
have  been  sAigled  out  for  special  treatment. 

The  functional  character  of  AV  is  determined  by  the  gravitational  field,  and  by  aero¬ 
dynamic  lift  and  drag  when  the  vehicle  is  within  the  atmosphere.  To  be  perfectly  general, 
we  may  include  continuous  thrust  as  a  limiting  condition  obtained  by  allowing  n  to  approach 
infinity. 

The  complete  set  of  AV's  corresponding  to  ell  possible  values  of  the  independent  state 
variables  comprises  a  (4n  -  6)-fold  infinity  of  solutions.  But  for  the  various  missions 
we  are  in  general  interested  only  in  those  satisfying  the  (4n  -  8)  equations, 


h  (AV)\ 

\apiJ 


const.  tL, 


(n-2) 


i.e, ,  only  those  solutions  for  which  AV  has  been  minimized  with  respect  to  all  the  state 
variables  except  tL  and  tft.  These  last  two  require  preferential  treatment;  in  fact,  one 
is  rarely  willing  to  allow  them  to  become  arbitrarily  large  in  the  process  of  minimizing 
AV. 


2.  The  Exact  Solution.  The  exact  solution  of  the  (4n  -  8)  equations  (II-2),  plus  the  aero- 
dynamic  and  Kepler ian  equations  of  motion  and  constraints,  is  an  extremely  difficult 
mathematical  task  even  with  unlimited  computational  facility,  much  less  for  a  self- 
contained  system.  It  has  been  carried  out  only  for  very  special  and  restricted  cases. 
However,  before  introducing  any  compromises  or  approximations,  it  would  be  well  to 
carry  the  logical  process  somewhat  further,  in  an  effort  to  predict  what  the  solutions  may 
actually  look  like,  and  what  one  would  do  with  them  if  they  were  available.  This  not  only 
provides  some  insight  into  the  design  of  an  operational  system,  but  it  may  also  suggest 
what  the  nature  of  the  approximations  should  be. 

In  order  to  provide  a  graphic  model,  we  will  consider  the  case  of  launch  from  earth  to  an 
arbitrary  orbital  rendezvous,  although  the  principles  deduced  may  be  generalized 
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further.  The  solution  of  equations  (II-2)  comprises  a  two-fold  infinity  of  solutions, 
corresponding  to  all  possible  values  of  t l  and  t^  Each  of  these  solutions  and  its 
corresponding  AV  may  be  represented  as  a  point  in  the  (tL,  t^)  plane.  Points  of  equal 
AV  may  be  joined,  forming  uniform  AV  contours.  The  resulting  contour  map  could  look 
like  that  of  Figure  2.  All  points  in  this  plane  represent  minimum  AV  with  respect  to  the 
state  variables  Pj;  but  tL  and  t^  are  retained  as  independent  variables.  Presumably  there 
would  be  at  least  one  ’’pit”  in  the  contour  relief  map,  corresponding  to  that  combination  of 
tL  and  t^  for  which  AV  approaches  a  to^al  minimum.  The  sense  of  the  ’’fall  line”  in  this 
figure  indicates  decreasing  AV. 

The  (4n  -  6)-fold  infinity  of  solutions  has  been  reduced  to  a  two -fold  infinity.  However, 
still  further  reduction  is  possible.  Except  for  certain  considerations  which  will  be  dis¬ 
cussed  subsequently.  Launch  time  need  not  in  general  be  preserved  as  a  free  parameter. 
What  the  operations  officer  is  really  concerned  about  in  terms  of  most  favored  solutions 
is  time  of  rendezvous  and  AV  required.  In  principle  he  does  not  care  when  he  has  to 
launch,  as  long  as  it  is  not  in  the  past.  Thus  the  following  additional  minimization  con¬ 
ditions  may  be  imposed. 


a  (AV) 


a  t. 


=  0 


const,  t 


R 


(n-3) 


The  set  of  solutions  satisfying  this  condition  is  represented  by  the  curve  AB  in  Figure  3. 

It  is  obtained  by  drawing  horizontal  lines  through  all  possible  values  of  t^,  and  selecting 
points  of  tangency  with  the  smallest  AV  contour.  As  one  moves  backward  in  time,  the 
curve  is  interrupted  by  the  vertical  line  BC  correjponding  to  launch  in  the  present.  Thus 
the  complete  set  of  solutions  availalle  in  trading  rendezvous  time  for  AV  is  the  curve 
ABC.  (There  is  no  merit  in  continuing  this  same  curve  beyond  A.  since  this  would  entail 
increasing  both  AV  and  t^,  which  is  a  bad  bargain.  After  "the  present”  reaches  point  A, 
the  desired  solutions  fall  on  the  vertical  "present”  line.) 

If  the  mission  prescribes  a  fuel  budget,  the  intersection  of  the  curve  ABC  with  the 
appropriate  AV  contour  is  the  required  solution-  If  one  is  prepared  to  sacrifice  fuel  to 
accomplish  earlier  rendezvous,  the  trade -off  is  effected  by  moving  down  the  curve  as  far 
as  fuel  considerations  permit.  The  only  trajectories  which  are  to  be  considered  fall  in 
this  singly  infinite  set,  and  no  other  points  in  the  tjj)  plane  need  be  computed. 

It  should  be  noted  that  the  curve  AB  may  be  obtained  alternately  by  the  following  condition: 


a  t 
Ft 


R 

L 


const.  AV 


0 


(n-4) 


This  demonstrates  that,  in  spite  of  the  fact  that  the  solutions  were  obtained  by  minimizing 
AV  we  have  simultaneously  solved  the  minimum  rendezvous  time  problem.  Thus  AV  and 
tR  do  not  prescribe  two  entirely  different  mathematical  approaches,  corresponding  to 
different  mission  requirements.  If  the  equations  have  been  properly  formulated,  one 
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TIME  OF  LIFT-OFF 


Figure  2.  Della  V  Contours  in  the  Rendezvous  Time, 
Liftoff  Time  Plane 


'  TIME  OF  LIFT-OFF 

THE  PRESENT 


Figure  3.  Rendezvous  Time  -  Trade-off  Curve 


common  computational  procedure  may  be  used  for  many  different  applications.  Only  in  the 
final  stage  of  making  a  logical  decision  is  it  necessary  to  introduce  the  criteria  appropriate 
to  the  particular  mission. 

The  contours  of  Figure  3  are  conceptual.  No  attempt  has  been  made  above  to  predict  their 
complete  representation.  It  would  therefore  be  useful  at  this  point  to  become  romewhat 
more  specific. 

The  first  observation  that  can  be  made  is  that  half  the  (tL,  t^)  plane  is  a  forbidden  region. 
There  is  a  45  degree  barrier  through  the  origin,  corresponding  to  the  fact  that  rendezvous 
can  never  precede  launch.  In  fact,  as  we  approach  this  barrier  the  AV  contour  lines  be¬ 
come  bunched  together;  a  relief  model  in  which  the  height  of  the  surface  represents  AV 
would  approach  infinite  steepness  along  this  45  degree  wall,  corresponding  to  infinite  AV 
for  rendezvous  in  zero  time. 

Another  observation  to  be  made  is  the  periodic  character  of  the  minimum  AV  regions, 
corresponding  to  favorable  launch  and  rendezvous  orientations  which  tend  approximately 
to  repeat  themselves  in  time.  There  will  be  two  such  periodicities  in  the  problem,  one 
corresponding  to  the  motion  of  the  launch  site,  and  the  other  that  of  the  rendezvous  target. 

The  fact  that  there  is  a  daily  launch  time  which  is  favorable  with  respect  to  AV  is  fairly 
obvious,  since  this  is  the  time  of  closest  approach  (or  crossover)  between  the  target 
orbital  plane  and  the  launch  site.  It  is  somewhat  less  obvious  that  there  is  likewise  a 
favorable  rendezvous  time.  However,  it  can  be  shown  analytically,  and  has  been 
demonstrated  by  simulation  studies  (see  Reference  9),  that  the  most  favorable  (i.e. , 
minimum  AV)  point  of  entry  into  the  target  orbit  tends  to  be  near  the  line  of  nodes,  which 
is  the  intersection  of  the  parking  and  target  orbital  planes.  Since  the  target  crosses  this 
line  twice  per  period,  it  introduced  AV  periodicity  with  respect  to  t^. 

In  the  case  of  launch  into  low  earth  orbit,  the  launch  site  position  will  have  a  daily  period 
corresponding  to  earth’s  rotation,  while  the  target  trajectory  will  have  a  much  shorter 
orbital  period.  The  effect  on  AV  is  depicted  in  Figure  4.  The  contour  ’’pit”  is  now  seen 
to  recur  periodically  with  respect  to  t^  as  the  target  site  rotates  every  couple  of  hours. 

The  daily  earth’s  rotation,  on  the  other  hand,  manifests  itself  as  a  slow  periodicity  with 
respect  to  tL,  introducing  long  horizontal  valleys  in  the  contour  map.  These  tend  to  re¬ 
peat  themselves  every  day,  but  Figure  4  follows  this  for  only  about  half  of  a  daily  period 

in  tL' 

As  the  AV  valleys  encroach  upon  th**  45  degree  causality  barrier  they  become  foreshortened, 
terminating  in  steep  cliffs.  Thus  the  daily  periodicity  is  interrupted  by  this  barrier.  If 
the  diagram  of  Figure  4  were  extended  up  into  higher  values  of  t^,  the  diurnal  periodicity 
with  respect  tc  tL  would  manifest  itself  before  being  terminated  by  the  barrier. 

Even  though  we  have  not  here  solved  the  exact  problem,  or  derive  precisely  the  shapes 
of  the  contours  for  the  multiburn  problem,  it  is  thus  possible  to  anticipate  certain 
symmetries,  periodicities,  and  limits  in  the  actual  solution.  The  trade-off  curves, 
designated  as  AB  in  the  previous  figure,  and  containing  the  family  of  most  favorable  solu¬ 
tions,  is  now  repeated  for  each  valley.  If  we  can  compute  solutions  along  the  repeated 
solid  line  segments  AB,  we  have  enough  information  for  making  the  correct  decision,  as 
well  as  the  appropriate  flight  plan. 
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TIME  OF  RENDEZVOUS 


Figure  4.  Contours  and  the  Approximate  Solution  for  a 
Typical  Rendezvous  Problem 
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Thus  we  have  demonstrated  that  a  common  mathematical  formulation  may  in  principle  be 
employed  for  different  problems,  missions,  and  objectives;  in  generating  a  single  family 
of  solutions  from  the  multi-dimensional  hyperspace  of  all  possible  trajectories,  we  have 
reduced  the  requirement  of  mission  flexibility  to  that  of  making  logical  decisions.  The 
mathematics  of  the  problem  calls  for  the  minimization  of  AV,  and  the  favorable  curves 
AB  may  be  generated  by  considering  various  values  of  t^,  and  minimizing  AV  with  respect 
to  all  other  variables. 

3.  The  Approximate  Solution.  There  remains  of  course  the  matter  of  obtaining  a  useful 
and  reliable  solution  to  this  mathematical  problem.  To  do  this  within  the  confines  of  a  self- 
contained  system  and  reasonably  finite  time,  it  is  clearly  necessary  at  the  present  time  to 
introduce  approximations.  However,  since  we  have  some  idea  of  the  form  of  the  solution, 
it  is  possible  to  make  approximations  and  simplifications  judiciously. 

First  we  note  that  for  the  low  earth  orbits  represented  in  Figure  4,  savings  in  Tr  obtained 
along  one  of  the  AB  segments  are  at  most  only  of  the  order  of  one  half  orbit.  Although  it 
is  conceivable  that  certain  situations  may  eventually  be  encountered  in  wh'f  .1  one  is  willing 
to  expend  extra  fuel  in  order  to  gain  only  a  few  minutes  of  rendezvous  time,  it  is  more 
likely  that  the  useful  trade-off  will  involve  hours  or  days,  rather  than  minutes.  In  other 
words,  if  there  is  to  be  any  significant  time  -  AV  trade-off,  it  will  be  not  so  much  where 
along  a  particular  AB  curve  to  find  the  desired  solution  as  which  of  the  AB  curves  to  select. 

Furthermore,  restriction  to  flight  plans  falling  precisely  on  the  AB  segments  results  in 
being  constrained  to  narrow  launch  windows,  since  there  are  large  gaps  in  usable  tL, 
corresponding  to  regions  of  excessive  AV.  Hold  problems  in  an  actual  launch  process, 
as  well  as  various  possible  military  requirements,  tend  to  make  wide  launch  windows  an 
important  consideration.  An  example  of  this  is  the  failure  of  the  Gemini  9  mission  of 
June  1,  1966,  which  missed  its  six-minute  launch  window  because  of  a  hold  at  T  minus 
three  minutes  resulting  in  a  two -day  delay. 

These  factors  lead  us  rather  to  look  for  solutions  following  along  the  AV  valleys,  shown 
as  a  series  of  dotted  steps  in  Figure  4.  This  series  of  solutions  is  obtained  by  considering 
various  values  of  tL,  and  minimizing  AV  with  respect  to  tR,  i,  e. , 


a  (AV)' 

lR 


const,  t^ 


0 


(D-5) 


The  jump  to  the  subsequent  step  or  "valley"  is  obtained  as  tL  increases  whenever  a  AV  is 
x'eached  which  exceeds  the  fuel  budget  for  a  particular  booster -vehicle  combination.  This 
series  of  steps  now  represents  our  trade-off  between  rendezvous  time  and  AV. 

It  is  these  AV  valleys  and  associated  wide  launch  windows  which  are  the  computational 
objectives  ot  the  Variable  Point  Guidance  and  Targeting  Program  (described  briefly  in 
Paragraph  5  of  this  section).  The  targeting  (or  flight  plan  generating)  portion  of  the  pro¬ 
gram  varies  tL  and  implicitly  seeks  to  minimize  AV  with  respect  to  the  other  state 
variables. 

It  should  be  noted  that  while  the  preceding  discussion  has  suggested  that  the  problem  is 
being  solved  from  earth  launch,  actually  *l  may  be  regarded  as  starting  from  any  initial 
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trajectory.  In  fact,  the  computation  is  constantly  being  repealed  (two  seconds  was  used  as 
the  cycle  time  in  simulation)  so  that  the  process  represents  not  merely  a  targeting  program, 
but  a  guidance  scheme  as  well.  In  other  words,  the  fact  that  at  any  instant  of  computation 
the  vehicle  may  be  halfway  to  its  target  merely  requires  the  inhibition  of  certain  logical 
steps,  such  as  those  involving  the  actual  launch  process,  range  safety  constraints,  etc.; 
the  basic  mathematics  remains  essentially  unchanged. 

The  exact  solution  of  the  basic  mathematics  would  require  solving  simultaneously  not  only 
ballistic  and  aerodynamic  equations,  but  also  all  the  minimization  and  constraint  require¬ 
ments  as  well.  The  most  obvious  simplification  is  to  break  this  complex  problem  into 
several  lesser  ones,  and  find  solutions  for  each  of  the  subsidiary  problems.  Clearly  the 
resulting  pieced-together  solutions  will  not  be  identical  to  the  overall  optimum,  but  if  the 
subdivision  is  accomplished  properly,  there  should  be  enough  uncoupling  to  make  the  re¬ 
sults  a  good  approximation. 

The  first  step  is  to  uncouple  the  pitch  and  yaw  guidance  through  the  atmosphere  as  a 
separate  entity,  with  the  exception  of  the  determination  of  launch  azimuth.  In  other  words, 
a  programmed  vertical  rise,  roll  to  proper  azimuth  and  pitch  maneuver  to  the  end  of 
sensible  atmosphere  is  employed.  Insertion  into  a  low  parking  orbit  is  accomplished  by 
the  generated  guidance  commands.  The  azimuthal  direction  of  the  parking  orbital  plane 
is  so  chosen  as  to  minimize  its  dihedral  angle  with  the  target  orbital  plan  (subject,  of 
course,  to  range  safety  constraints,  etc.).  This  is  the  extent  to  which  the  launch  process 
is  coupled  to  the  overall  requirements. 

4.  Guidance 


a.  Atmospheric  Exit  --  Open  Loop  Guidance.  During  the  initial  portions  of  the  launch 
trajectory  the  constraints  imposed  by  aerodynamic  forces  are  so  severe  as  to  outweigh 
consideration  of  optimizing  propellant  utilization.  While  the  launch  vehicle  is  in  the  lower 
portion  of  the  atmosphere,  therefore,  the  actual  trajectory  cannot  depart  greatly  from  a 
gravity  turn  without  exceeding  the  structural  limits  of  current  booster  configurations. 

Since  the  bending  moment  on  the  vehicle  stack  is  essentially  proportional  to  the  product 
of  angle  of  attack  and  dynamic  pressure,  the  angle  of  attack  must  be  kept  small,  particularly 
over  the  duration  of  maximum  dynamic  pressure.  The  dynamics  of  stage  separation  also 
require  a  small  angle  of  attack  if  appreciable  dynamic  pressure  remains  at  this  time. 

These  effects  are  further  exaggerated  if  the  final-stage  vehicle  is  itself  a  lifting  body. 

Thus,  guidance  during  this  phase  is  generally  considered  to  be  best  handled  by  a  program 
which  will  achieve  trajectories  close  to  gravity  ;’urns.  There  are  two  parameters  which 
can  be  varied  (within  the  limits  imposed  by  range  safety  to  optimize  the  trajectory  during 
the  atmospheric  phase.  They  are  launch  azimuth  and  pitch-over  gain  (or  other  equivalent 
gravity  turn  parameter).  Once  the  overall  flight  planning  routine  has  selected  these 
parameters  properly,  and  included  an  optimum  trajectory  beyond  the  atmosphere,  further 
optimization  of  the  trajectory  within  the  severe  angle  of  attack  limitations  will  yield  little 
improvement  in  overall  efficiency. 

The  attitude  control  during  this  phase  must  operate  to  satisfy  the  angle  of  attack  limits  in 
spite  of  gust  loads  and  wind  shear  effects.  This  will  usually  eliminate  the  possibility  of 
flying  the  nominal  trajectory  with  tight  loop  attitude  control  and  requires  aerodynamic 
stability  augmentation.  The  nature  of  the  control  loop  depends  directly  on  the  vehicle 
characteristics,  including  aerodynamic  parameters,  structural  am1  heating  limits,  bending 
mode  shapes,  sensor  locations  and  center  of  gravity  motion  during  burning  (including  the 
effects  of  propellant  slosh  in  liquid  propellant  stages). 
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b.  Closed-Loop  Guidance.  From  the  point  of  view  of  overall  guidance  the  resultant 
effect  of  atmo spker ic  disturbanc e s  and  off  nominal  vehicle  performance  will  be  a  deviation 
of  the  vehicle  stage  from  the  nominal  upon  exit  from  the  effective  atmosphere.  During  the 
subsequent  portion  of  the  launch  phase  the  desired  cutoff  conditions  must  be  achieved  in  a 
near  optimum  manner.  This  is  the  function  of  the  closed-loop  guidance  used  throughout  a 
burn  phase.  The  ascent  phase  guidance  is  greatly  simplified  if  all  missions  use  a  low 
altitude  circular  orbit,  either  as  the  primary  orbit  for  the  mission,  or  as  a  parking  orbit. 
In  this  case,  a  restricted  range  of  standard  ascents  is  derived  for  all  missions  with  a  given 
booster,  the  relatively  small  variations  with  launch  azimuth,  payload  weight  and  orbit 
altitude  being  dealt  with  in  the  latter  portion  of  the  ascent  phase. 


The  simplest  method  of  guidance  for  this  phase  which  is  not  restricted  to  a  limited  range 
of  orbits  specifies  the  velocity  vector  as  a  function  of  position  at  cutoff.  This  method  is 
used  for  two-burn  entry  into  orbit,  in  the  current  Variable  Point  Guidance  and  Targeting 
Program. 

In  this  case,  a  "velocity  required"  is  defined  at  each  instant  of  flight  and  the  vehicle  thrust 
vector  oriented  to  reduce  the  velocity-to-be -gained.  Pitch  and  yaw  steering  is  used  to 
maintain  the  velocity  to  be  gained  vector  non -rotating  in  inertial  space.  Cut-off  occurs 
when  velocity-to-be-gained  becomes  a  null  value.  Experience  has  shown  that  his  is  an 
efficient,  though  not  strictly  optimal  method  of  guidance. 

5.  Orbital  Transfer  and  Rendezvous.  The  problem  of  orbit  transfer  and  rendezvous  is  to 
make  a  minimum  AV  transfer  from  an  initial  orbit  (usually  a  low  parking  orbit,  not 
necessarily  circular)  to  a  specified  target  orbit,  and  to  rendezvous  in  phase. 

The  parking  orbit,  from  which  this  phase  of  the  problem  begins,  has  already  been  optimized 
with  respect  to  azimuthal  direction,  to  the  extent  that  the  plane  of  this  orbit  makes  a 
minimal  dihedral  angle  with  the  target  orbital  plane;  this  reduces  the  amount  of  plane  change. 
Strictly  speaking,  this  is  only  approximately  a  valid  criterion  for  minimizing  AV.  However, 
simulation  studies  show  it  to  be  an  excellent  approximation,  and  it  accomplishes  the  avowed 
purpose  of  uncoupling  the  launch  process  from  the  rest  of  the  pr  blem. 

The  first  approximation  to  make  in  solving  the  problem  of  orbital  transfer  and  rendezvous 
from  a  parking  orbit  is  to  restrict  the  number  of  impulse  burns.  (The  impulse  is  of  course 
itself  an  approxi  .ation  which  may  be  shown  to  be  very  good  for  all  burns  except  launch.) 

Here  we  set  n  =  3.  However,  the  computer  solves  this  3 -burn  problem  by  looking  for  the 
minimum  AV  2 -burn  solution  necessary  to  make  an  orbital  transfer  without  rendezvous. 

This  same  AV  is  then  imposed  on  a  3 -burn  solution  which  actually  effects  rendezvous. 


The  minimum  AV  2-burn  trajectory  to  effect  rendezvous  would  involve  4n  -  6  =2  free 
variables.  However,  since  here  we  are  solving  the  problem  as  though  it  were  orbital 
transfer  without  rendezvous,  the  target  pjjase  is  likewise  left  free,  so  that  altogether  there 
are  three  free  variables.  The  desired  AV  min  solution  with  respect  to  these  3  variables 
causes  the  vehicle  to  arrive  at  the  target  orbit  with  a  velocity  7enter>  out  of  phase  with  the 
target.  Instead  of  giving  the  vehicle  the  velocity  the  target  has  at  this  point,  we  insert  with 
a  velocity  7phase  which  causes  the  vehicle  to  overtake  the  target  after  one  or  more  orbital 
periods  and  arrive  back  at  this  point  precisely  in  phase,  whereupon  the  third  and  last  burn 
gives  it  the  target  velocity 


Since: 


7  =7  -7  =(7  -7  )  +  (7  -7  ) 

min  target  enter  phase  enter  target  phase 
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We  can  decompose  the  second  burn  of  the  minimum  AV  solution  for  2 -burn  orbital  transfer 
into  two  AV  impulses,  the  first  of  which  accomplishes  phasing  and  the  second  rendezvous, 
but  with  the  same  total  AVmjn.  Thus,  orbital  transfer  is  effectively  accomplished  with 
minimum  2-burn  AV,  and  phasing  into  rendezvous  with  zero  AV.  This  procedure  makes 
the  fine  phasing  adjustments.  Larger  corrections  are  accomplished  more  rapidly  while  in 
low  parking  orbit. 

The  program  for  solving  the  minimum  AV  2 -burn  transfer  problem  by  searching  over  three 
state  variables  has  been  programmed  and  tested  on  various  target  orbits.  However,  it  has 
been  observed  that  the  optimum  point  for  transferring  to  target  orbit  always  falls  very 
close  to  the  line  of  nodes,  which  is  the  intersection  of  the  parking  and  target  orbital  planes. 
If  we  introduce  this  as  an  arbitrary  constraint,  and  add  further  the  requirement  that 
transfers  from  parking  orbit  to  transfer  orbit,  and  transfer  orbit  to  target  orbit,  are  always 
to  be  tangent  transfers,  then  we  have  eliminated  all  three  free  variables.  Previous  studies 
have  shown  the  results  of  this  simplifying  approximr 'ion  to  be  very  close  to  those  of  the 
exact  minimization  process;  hence,  the  simplified  ^  ;sion  has  been  incorporated  into  the 
Variable  Point  Guidance  and  Targeting  Program. 

The  solution  of  the  rendezvous  problem  described  here  has  involved  approximation  by  sub¬ 
division  into  phases;  namely,  launch,  orbital  transfer,  and  rendezvous.  Now  one  unifying 
final  step  is  introduced;  the  plane  changes  are  redistributed  over  the  nodal  burns  in  an 
optimal  manner.  The  computation  minimizes  AV  with  respect  to  each  of  the  plane  changes, 
subject  to  the  constraint  that  the  total  plane  change  is  conserved.  This  adjustment  makes 
the  last  correction  in  the  approximate  solution  of  the  minimum  AV  problem. 

After  lift-off,  when  the  engines  are  on,  the  guidance  commands  are  derived  from  the  flight 
plan  requirements.  After  each  burn  phase  and  as  the  mission  progresses  the  consideration 
of  various  computational  alternatives  is  suppressed,  in  accordance  with  the  fact  that  the 
existing  state  of  the  system  is  the  result  of  previous  logical  commitments.  Terminal 
guidance  or  the  closure  mode  introduces  no  changes  in  the  guidance  equations  themselves, 
except  for  the  time  or  distance  between  burns  and  a  coordinate  transformation  which  con¬ 
verts  on-board  sensor  data  from  relative  coordinates  to  inertial  coordinates.  Mid-course 
and  terminal  maneuvers  are  thus  reduced  to  a  special  case  of  orbital  transfer  and  rendez¬ 
vous. 
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SECTION  HI 


EQUATIONS  AND  LOGIC 


1.  Introduction.  The  purpose  of  this  station  is  to  describe  how  the  Variable  Point  Guid¬ 
ance  and  Targeting  technique  works.  The  basic  equations  and  logic  as  well  as  the  required 
inputs  and  outputs  of  the  system  are  given  and  discussed.  The  derivation  of  the  equations 
are  given  in  Appendix  I.  The  guidance  and  targeting  requirements  for  the  various  phases  of 
the  rendezvous  mission  from  before  launch  to  the  beginning  of  the  station  keeping  phase  are 
considered. 

The  phases  of  the  mission  are  given  in  Table  IV.  A  summary  cf  the  computational  opera¬ 
tions  performed  during  the  various  phases  of  a  typical  rendezvous  mission  are  also  included. 

rXhe  flight  plans  used  by  Variable  Point  Guidance  and  Targeting  are  selected  from  among 
several  possible  plans,  involving  the  use  cf  an  ascent  burn,  plus  any  or  all  of  the  following: 
parking  orbits,  transfer  orbits,  plane  change  maneuvers,  small  corrections  maneuvers, 
and  an  adaptation  maneuver  into  the  target  orbit.  Although  the  number  of  burns  used  is 
discrete  (n  -  2),  the  flight  plans  span  a  continuum  in  the  sense  that  they  tend  to  vary  con¬ 
tinuously  with  time.  The  type  of  flight  plan  to  be  used  in  a  given  situation  is  selected  in  the 
Variable  Point  Guidance  and  Targeting  computation,  based  on  the  instantaneous  geometry  of 
target  and  vehicle.  The  vehicle  can  be  guided  to  any  target  whose  orbit  is  within  the 
physical  capability  of  the  vehicle’s  propulsion.  Last-minute  changes  in  vehicle  configuration 
(payload  weight,  propellant  load,  etc. )  can  be  accepted  with  a  minimum  amount  of  repro¬ 
gramming  and  corrections  for  weight  and  thrust  deviations,  winds,  variation  in  lift-off  time, 
steering  lags,  etc. ,  can  be  made  in  flight. 

The  scheme  requires  a  knowledge  of  the  vehicle's  position  and  velocity,  which  can  be  ob¬ 
tained  from  data  supplied  by  any  form  of  inertial  measurement  unit  IMU  ,  inertial  platforms, 
strapped  down  inertial,  stellar-monitored  inertial,  etc. ,  are  all  compatible  with  Variable 
Point  Guidance  and  Targeting.  A  navigation  error  analysis,  which  necessarily  must  be 
performed  for  any  hardware  configuration,  depends  primarily  on  the  inertial  and  computing 
devices,  and  not  on  Variable  Point  Guidance,  whose  mathematical  errors  are  small  by 
comparison.  The  Variable  Point  Guidance  equations  can  be  solved  in  the  ground  based 
computer  of  a  radio-command  system,  provided  the  measurement  of  vehicle  position  and 
velocity,  and  the  communication  link  to  the  vehicle  can  be  maintained. 


Table  IV.  Phases  of  Rendezvous  Mission 


1.  Prelaunch 

Targeting  and  Flight  Planning 

2.  Lift-off  Through  Boost 
to  First  Cut-oif 

Boost  Guidance  Steering,  Staging, 
Navigation,  Cut-o'f  Computation 

3,  Coast  and  Parking  Orbit 

Orientation  Control,  Flight  plan  Up¬ 
dating,  Time  to  Start  Next  Burn, 
Navigation 

A.  Transfer  Burn 

Guidance  Steering,  Staging,  Naviga¬ 
tion,  Cut-off,  Computation 

5.  Coast 

See  (3) 

6.  Phasing  Control  Burn 

See  (4) 

7.  Phasing  Orbit 

Continue  (3)  until  Closure  Mode 
(Terminal  Phase*) 

Closure  Mode  (Terminal) 

Series  of  Burn- Coast  Phases  to 
Rendezvous  Point 

Burn 

Navigation  and  Guidance  Steering 
and  Cut-off 

Coast 

Navigation  and  Orientation  Control 
with  Start  Burn  Computation  for 
Shorter  Segments 

♦Terminal  guidance  usually  refers  to  the  phase  of  a  mission  in  which 
onboard  sensors  (such  as  radar)  are  used  to  acquire  a  target  and 
furnish  information  which  can  be  processed  and  used  to  provide 
guidance  signals  for  further  corrective  action. 


The  prime  functions  of  the  Variable  Point  Guidance  and  Targeting  computation  in  the  pre¬ 
launch  phase  are  navigation  (position  and  velocity),  targeting  and  flight  planning  and  pro¬ 
viding  information  to  the  launch  conductor.  The  targeting  function  is  performed  by  genera¬ 
ting  the  flight  plan  from  which  data  can  be  fed  to  a  display  for  the  launch  conductor.  In  the 
present  Variable  Point  Guidance  and  Targeting  formulation,  the  launch  conductor  commands 
liftoff. 

After  liftoff,  the  prime  functions  are  navigation,  targeting,  guidance,  steering,  and  issuance 
of  discretes.  The  navigation  and  targeting  functions  are  performed  continuously.  Guidance 
determines  the  present  action  needed  to  implement  the  flight  plan  and  translates  this  into 
vehicle  commands.  The  scheme  is  based  on  a  restart  capability  in  stages  beyond  those  that 
are  burned  to  depletion  in  the  first  burn.  This  is  essential  to  gain  the  fullest  and  most 
efficient  use  of  the  propulsion  capability  of  the  vehicle  to  rendezvous  with  targets  in  the 
widest  variety  of  orbits.  It  is  primarily  a  mission  requirement  and  not  a  guidance  limitation. 
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The  inputs  to  Variable  Point  Guidance  and  Targeting  are  shown  in  Figure  5.  The  target 
ephemeris  is  the  only  target  data  which  need  be  furnished  to  the  computer  prior  to  the 
terminal  phase.  Typical  IMU  inputs  are  increments  of  velocity  and  gimbal  angle.  A  "Go- 
inertial"  signal,  signifying  that  the  IMU  has  entered  its  operating  mode,  causes  the  normal 
navigation  computation  to  begin.  The  time  synchronization  signal  references  the  computer 
clock  to  Greenwich  mean  time.  The  lift-off  signal  instructs  the  computer  to  take  ever  active 
control  of  the  vehicle.  The  vehicle’s  AV  budget  is  supplied  for  the  purpose  of  determining 
launch  windows.  An  instant  of  time  is  within  a  launch  window  if  the  AV  required  for  rendez¬ 
vous  (launching  at  that  instant)  is  less  than  the  AV  budget.  A  pitch  program  and  multiplier 
suitably  defined  for  the  particular  vehicle  configuration  are  furnished  to  assure  compliance 
with  the  heating  and  structural  constraints  in  the  atmospheric  exit  phase. 

The  outputs  furnished  to  the  launch  conductor  during  the  prclaunch  phase  are  indicated  in 
Figure  5.  The  outputs  after  lift-off  include  steering  commands  during  both  thrust  and  coast 
periods;  engine  ignition  and  cutoff;  and  sequencing  discretes  for  staging,  search,  etc.  The 
steering  commands  are  typically  angular  rate  commands  about  each  of  the  three  body  axes. 
Precision  engine  cut-off  commands  are  issued  at  the  end  of  each  burn.  Ignition  commands 
are  given  to  initiate  burns  or  to  complete  burns  as  required  following  the  jettisoning  of  a 
depleted  stage.  Staging  commands  are  issued  when  the  propellant  of  a  givr.A  stage  is 
depleted.  Search  is  initiated  when  the  predicted  relative  range  to  the  target  is  within  the 
terminal  sensor  capability  Flight  plans  are  generated,  each  cycle  before  lift  off  and 
revised  periodically  during  coast  periods. 

Guidance  commands  are  generated  each  cycle  during  a  burn  phase  based  upon  a  generally 
fixed  flight  plan.  Major  changes  in  the  target  orbit  or  changes  in  the  mission  can  be  in¬ 
serted  any  time,  however,  in  the  present  program  the  modified  flight  plan  for  a  major 
change  is  initiated  only  during  a  coast  period.  Minor  changes  such  as  ephemeris  updating 
can  be  accommodated  at  any  time.  In  the  terminal  phase  where  radar  data  would  be  used  to 
refine  the  target  ephemeris,  the  flight  plan  is  modified  between  each  burn. 

2.  Targeting  and  Flight  Planning.  The  types  of  flight  plans  established  by  the  present 
Variable  Point  Guidance  and  Targeting  logic  handle  the  complete  three-dimensional  rendez¬ 
vous  problem,  in  which: 

1.  )  The  launch  site  generally  lies  outside  the  target  orbit  plane 

2.  )  The  target  is  in  any  arbitrary  elliptical  earth-orbit  with  any  phase  in 

its  orbit. 

3.  )  Launch  can  occur  at  any  time. 

Flight  plans  are  built  up  from  segments  involving  ascent  into  a  low-altitude  parking  orbit, 
transfer  orbit,  phasing  orbit,  plane  change  maneuver,  and  entry  into  the  target  orbit  with 
corresponding  equivalent  position  and  velocity. 

The  use  of  two  burn  direct  ascent,  in  which  the  time  constraint  is  satisfied  by  trajectory 
phasing  is  not  at  present  included  since  the  launch  window  for  this  method  is  relatively 
narrow.  The  direct  intercept  case  falls  out  as  a  singular  point  in  the  present  program-  If 
time  of  lift-off  is  controlled  closely  enough  a  direct  ascent  solution  will  result.  In  the 
present  program  pulsing  orbits  and  full  orbit  phasing  are  used  as  the  primary  flight  plan; 
bielliptic  flight  plans  are  used  only  in  the  closure  mode,  that  is  only  during  the  last  360 
degree  of  flight  before  the  planned  rendezvous  point. 
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Figure  5.  Variable  Point  Guidance  and  Targeting  Block  Diagram 


22 


The  mission  starts  with  ascent  into  a  nominally  circular  low-altitude  parking  orbit,  the 
second  burn  being  90  degree  downrange  from  the  launch  site.  The  vehicle  spends  0,  1/2, 

1,  1  1/2  . . .  periods  in  the  parking  orbit  until  the  desired  relative  phase  between  target  and 
vehicle  is  grossly  established.  Then  a  third  burn,  initiates  a  transfer  to  the  target  orbit 
altitude,  reaching  the  target  orbit  at  the  line  of  nodes,  as  shown  in  Figure  6.  The  shape  of 
the  transfer  orbit  is  chosen  so  that,  if  the  plane  of  the  parking  orbit  were  rotated  about  the 
line  of  nodes  into  the  plane  of  the  target  orbit,  the  transfer  would  be  tangent  at  the  initial 
and  final  points.  Since  both  nodes  are  considered  as  potential  rendezvous  point  this  policy 
gives  close  to  the  optimum  transfer  for  one  of  the  nodes  even  for  final  orbits  of  relatively 
high  eccentricity  and  plane  change.  The  policy  generates  an  essentially  optimal  transfer 
when  the  final  orbit  line  of  apsides  is  coincident  v/ith  the  line  of  nodes.  A  fourth  burn  at  the 
target  altitude  initiates  a  full  phasing  orbit  whose  major  axis  is  adjusted  to  give  precise 
control  over  time  of  return  to  the  point  of  initiation.  A  fifth  burn  completes  the  rendezvous. 

Both  chase  and  lob  solutions  are  used  in  this  type  of  flight  plan.  In  chase,  the  phasing  orbit 
is  below  the  target  orbit,  and  in  lob,  it  is  above.  The  lob  type  ilight  plan  has  been  found  to 
be  useful  in  maintaining  wide  launch  windows  for  low  altitude  target  oruit  situations.  Low 
altitude  orbits  are  considered  to  be  those  below  about  1000  nautical  miles.  For  target  above 
1000  nautical  miles  phasing  can  be  effectively  accomplished  for  any  time  of  lift-off  below  the 
target  orbit.  For  targets  at  very  high  altitudes  (such  as  synchronous)  several  alternative 
phasing  solutions  may  exist  in  which  the  phasing  can  be  accomplished  below  the  target  orbit 
and  the  alternatives  increase  with  increasing  altitude.  It  has  been  observed  that  the  best 
solution  for  all  cases  is  to  maintain  phasing  as  near  to  the  target  orbit  as  possible  vdthout 
sacrificing  flight  time.  This  is  accomplished  by  proper  selection  of  the  number  of  parking 
orbits.  If  the  target  is  high  enough,  the  phasing  orbit  can  be  made  equal,  within  a  small 
time  increment  of  the  target  orbit  and  the  phasing  orbit  become  unnecessary.  These  are 
regions  about  the  singular  points.  The  synchronous  satellite  orbit  allows  fifteen  singular 
solutions.  The  lunar  orbit  allows  about  268  singular  solutions.  The  distinction  between 
high  and  low  orbits  is  made  in  the  program  by  testing  the  difference  between  the  target 
orbit  period  and  the  sum  of  the  transfer  orbit  period  and  the  parking  orbit  period. 

The  distinction  is  necessary  because  not  all  phasing  orbits  can  be  used  due  to  the  finite 
size  of  earth.  The  minimum  perigee  for  a  phasing  orbit  must  be  greater  than  or  equal  to 
the  radius  of  the  earth  plus  atmosphere.  Since  tangent  transfers  are  used  from  a  low 
altitude  parking  orbit  the  minimum  size  of  the  phasing  orbit  must  be  about  equal  to  the 
transfer  orbit. 

For  the  higher  targets  orbit  the  period  of  the  phasing  orbit  does  not  need  to  approach  the 
minimum  and  when  the  difference  is  large  enough  (between  parking  orbit  and  target  orbit) 
the  minimum  period  is  arbitrarily  set  at  one  parking  orbit  period  smaller  than  the  target 
orbit  period.  When  the  sum  of  the  transfer  orbit  period  and  the  parking  orbit  period  is 
greater  than  the  target  orbit  period,  the  phasing  orbit  must  be  between  the  transfer  orbit 
and  the  target.  If  the  target  orbit  is  exceptionally  low  several  phasing  orbits  may  be 
required  to  remain  beneath  the  target  orbit.  For  these  cases  it  is  preferable  for  most  of 
the  launch  windows  to  allow  the  phasing  orbit  to  lob  or  pass  outside  the  target  orbit. 

In  either  case  the  number  of  parking  oi^its  are  first  computed  based  upon  the  use  of  one 
minimum  phasing  orbit  for  each  node.  The  number  is  obtained  by  truncating  to  full 
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NOTE.  NUMBERS  REPRESENT  MAJOR  BURNS 


Figure  6.  Full  Orbit  Phasing 
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increment  the  minimum  solution.  The  phasing  orbit  is  then  trimmed  upward  for  a  usable 
solution.  The  solution  to  be  used  after  lift-off  is  chosen  from  the  two  possibilities  (one  at 
each  node). 

The  Greenwich  time  of  rendezvous  Tr  can  be  expressed  in  terms  of  target  motion  from  an 
initial  reference  time  and  a  target  flight  time: 


tr  =  Tt  +  nt  Tt 


(hi-  i) 


where  T^  is  the  Greenwich  time  of  the  next  arrival  of  the  target  at  a  prescribed  rendezvous 
point,  is  the  period  of  the  target  orbit  and  m  is  the  number  of  full  target  orbits  before 
rendezvous.  The  equation  (HI- 1)  accounts  for  the  fact  that  several  full  orbits  after  first 
time  of  arrival  of  the  target  at  the  selected  rendezvous  point  may  be  required  to  allow  the 
satellite  interceptor  to  arrive  at  the  same  time.  The  time  of  rendezvous  may  similarly  be 
expressed  in  terms  of  the  interceptor  flight  plan  as  follows: 
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where  T[  is  the  Greenwich  time  of  next  arrival  of  the  interceptor  at  a  prescribed  point, 

(the  point  of  the  entry  into  the  parking  orbit,  for  example)  n.  is  the  number  of  parking 
orbits,  r.  is  the  period  of  the  parking  orbit. 

The  rendezvous  time  is  found  by  solving  these  equations  simultaneously  in  an  iterative  proce¬ 
dure.  In  this  process,  Tt  andTt  are  found  from  the  target  ephemeris.  T[  is  fixed  by  the 
altitude  of  the  parking  orbit,  T^  is  known  from  the  characteristic  flight  time  of  the  payload 
booster  configuration  to  reach  and  enter  the  parking  orbit.  The  unknowns  are  the  number 
of  interceptor  parking  orbits,  the  size  and  number  of  phasing  orbits,  and  the  number  of 
target  orbits.  The  method  for  solving  the  equation  is  as  follows: 

First  solve  for  the  above  equation  with  nt  =  n.  =  0 

Assume  first  that  one  phasing  orbit  is  used  and  set  its  period  at  its 
limiting  value,  then,  solve  for  the  only  remaining  unknown: 
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(in- 3) 


(n  =  1  on  first  trial). 

0 

The  limiting  value  tl  is  equal  to  7X  the  transfer  orbit  period  for  low  altitude  target  orbits 
or  it  is  equal  to  rt  -  r,  for  high  orbits.  A  test  is  made  on  (rt  -  -  tx)  to  determine 

which  limit  to  use.  If  the  quantity  is  greater  than  zero  the  high  limit  is  used,  if  the 
quantity  is  negative  the  low  limit  is  used.  If  nj  is  negative  nt  must  be  increased.  When 
n^  is  positive,  it  is  then  truncated  and  the  period  r  is  computed. 
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T,  = 
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TJ  ~  Wt  ~  Ti 
n0 


(ni-4) 


If  r  ^  is  greater  than  r4  then  n  ,  is  increased  by  one  unit. 

0  l  0 


Similar  solutions  are  made  for  both  of  the  intersections  of  the  line  of  nodes  of  the  target. 
The  solutions  are  compared  and  one  is  chosen  for  the  flight  plan  for  that  particular  case. 
For  lifting  off  from  a  ground  station  account  is  taken  of  the  placement  of  the  line  of  nodes 
ty  azimuth  control  and  the  off-flight  to  the  ncie  from  lift-off  is  accounted  for  in  the  total 
time  of  arrival  at  the  rendezvous  point.  A  flow  chart  of  this  process  is  shown  in  Figure  7. 
The  max  mum  number  of  target  orbits  before  rendezvous  is  limited  by  an  input  setting, 
nTmax*  The  launch  director  can  vary  this  setting  depending  on  his  willingness  to  trade 
fuel  for  time.  If  he  wishes  to  obtain  the  greatest  fuel  efficiency  he  must  set  in  a  large 
enough  number  to  assure  that  only  chase  solutions  will  be  obtained.  (nmax  =  3  was  large 
enough  for  a  1000  mile  altitude,  circular  orbit  case).  Setting  in  a  smaller  number  will 
result  in  generating  lob  solutions  when  they  are  required. 

With  the  complete  planar  trajectory  laid  out,  the  velocity  and  flight  path  angle  entering  and 
leaving  every  burn  is  known.  This  data  is  the  input  to  the  plane  change  split  routine,  whose 
output  then  completes  the  solution.  The  problem  is  solved  twice,  once  for  a  trail  rendez¬ 
vous  point  at  each  end  of  LN,  providing  two  options  with  different  rendezvous  times  and/or 
AV  requirements. 

The  total  required  plane  change  is  removed  at  burns  2,  4  and  5.  A  pure  plane  change 
maneuver  is  made  at  point  "b"  or  "c"  if  its  inclusion  decreases  the  total  mission  AV.  The 
amount  of  plane  changes  removed  at  each  of  tiies e  burns  is  determined  by  optimally  split¬ 
ting  the  total. 

If  one  wanted  to  rendezvous  by  the  most  fuel- efficient  means  on  every  flight,  one  would 
restrict  the  flight  plan  to  launch  when  the  launch  site  lies  in  the  target  orbit  plane,  since 
AV  requirements  are  minimum  when  the  required  plane  change  is  zero.  The  vehicle 
should  then  be  launched  into  circular  parking  orbit  at  the  lowest  altitude  commensurate 
with  atmospheric  heating  and  drag  constraints.  It  should  remain  in  that  orbit  until  the 
phase  is  grossly  correct,  then  enter  a  Hohmann  transfer  to  the  apogee  of  the  target  orbit, 
and  complete  fine  control  over  phasing  with  a  full  orbit  tangent  at  the  target  orbit  apogee. 
The  program  produces  this  flight  plan  for  in  plane  launch  since  for  this  special  case  the 
line  of  nodes  is  arbitrarily  set  along  the  major  axis  of  the  target  orbit. 

The  efficiency  of  the  full-orbit  tangent- tangent  method  compared  to  an  optimum  nun- rendez¬ 
vous  trajectory  into  the  same  final  orbit  from  the  same  launch  site  location  ranges  from 
good  to  excellent,  depending  on  the  orbit  eccentricity,  the  angle  between  the  line 
of  nodes  and  the  line  of  apsides  and  the  angles  between  the  planes.  As  the  angle  between 
LA  and  LN  decreases  from  90  degree  to  0  degree,  ana  /or  as  the  eccentricity  of  the  target 
orbit  varies  from  unity  to  zero,  the  tangent  transfer  becomes  essentially  optimal  over  a 
wide  range  of  out- of- plane  angles. 

The  full-orbit  method  of  phasing  has  the  desirable  characteristic  that,  for  a  given  inertial 
location  of  the  launch  site,  and  a  given  final  target  orbit  to  be  attained,  the  final  time  of 
rendezvous  has  no  effect  on  AV.  That  is,  for  mission  situations  in  which  only  the  phase 
of  target  is  varied,  the  resulting  trajectories  all  use  the  same  AV,  but  rendezvous  occurs 
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Figure  7.  Decision  Flow  Chart  Full  Orbit  Phasing 


at  different  times.  The  decoupling  between  time  of  rendezvous  and  AV  is  possible  because 
the  burn  at  the  high-altitude  *md  of  the  transfer  orbit  can  be  split  into  two  collinear  sub¬ 
burns,  performed  an  orbit  apart,  the  sum  of  which  is  veetoriaily  equal  to  the  original  single 
burn.  The  AV  allotted  to  each  sub-burn  controls  the  size  and  period  of  the  phasing  orbit, 
hence  the  time  of  arrival  of  the  interceptor  at  the  specified  rendezvous  point.  The  logic 
which  implements  this  phasing  method  solves  the  phasing  problem  first,  and  the  plane- 
change  problem  second. 

a.  On-Orbit  Tangent  Transfer.  The  transfer  trajectory  is  established  as  that  which  is 
tangent  to  the  target  orbit  at  the  intersection  of  the  target  orbit  with  LN,  and  tangent  to  the 
circular  parking  orbit. 

The  parking  orbit  velocity  vector  at  the  node  is: 
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The  tangent  of  flight  path  angle  at  the  node  and  the  radius  vector  at  the  nodes  are  computed 
for  both  the  target  orbit  and  the  parking  orbit.  The  tangent  transfer  (see  Figure  8)  is 
determined  by  finding  range  angle  from  the  node  on  the  target  orbit  to  the  tangent  point  on 
the  parking.  There  is  one  such  angle  for  each  node  defined  by  its  inverse  tangent  half  angle 
function. 
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The  trigonometric  identity  for  the  half  angle  gives  the  sine  and  cosine  of  the  angle. 
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The  sign  and  cosine  of  the  angle  is  computed  from 
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Figure  8.  Transfers  From  Parking  Orbit  to  Target  Orbit 
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The  time  at  node  is  obtained  as  follows: 
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b.  Optimal  Plane  Change  Split.  It  was  stated  above  that  the  rendezvous  problem  was 
solved  by  reducing  to  a  planar  case  first,  solving  for  the  shape  of  trajectory 
legs  to  satisfy  the  time  constraint,  and  then  secondly  determining  the  plane  for  each  of 
these  legs.  A  means  to  divide  the  total  plane  change  to  be  made  into  optimally- sized 
portions  is  as  follows: 
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AVT  “2  AVi  (Hi-H) 
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where  AVp  is  the  total  mission  AV  which  is  to  be  minimized.  AV|,  the  velocity  increment 
at  the  ith  impulse  point,  may  be  written  as 
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where  A.  and  B.  are  defined  as: 

Ai  =  4=  +  VfL  '  2ViE  ViL  IE  Sin^  iL 
Bi  =  2ViE  ViL  C0Sr  IE  C0Sr  iL 


and  where: 

velocity  entering  (leaving)  the  impulse  point 
flight  path  angle  entering  (leaving)  the  impulse  point 

Plane  change  angle  taken  at  that  impulse  point 
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Coefficients  Aj  and  are  constants  for  che  ith  burn  determined  fully  by  the  in-plane  por¬ 
tion  of  the  rendezvous  problem  solution.  The  total  plane  change  to  be  taken  (cv^  is  the 
dihedral  angle  between  the  interceptor's  initial  plane  and  the  target  orbit  plane.  It  can  be 
written  as  the  sum  of  the  plane  change  made  at  each  burn: 
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The  problem  is  stated  as  one  of  finding  that  set  of  aj's  which  minimizes  AV^  as  given  by 
equations  (III- 14)  subject  to  the  constraint  equation  (HI- 18).  The  solution  to  this  problem  is 
readily  obtained  by  the  method  of  Lagrange  multipliers. 

A  necessary  condition  for  an  extremum  is  that  the  incremental  AV  for  an  increment  of 
plane  change  angle,  8  (AV)  ,  be  a  constant  be  a  constant  for  each  burn.  That  is: 
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Screening  techniques  have  been  developed  to  set  up  valid  arrays  of  aim  point  conditions  to 
be  used  by  the  plane  change  split  routine.  The  objectives  of  screening  are: 

1.  )  To  exclude  points  from  plane  change  processing  which  would  cause 

convergence  of  the  routine  to  a  maximum  rather  than  a  minimum 

2.  )  To  cut  down  computation  time  caused  by  slowly  converging  solutions 

3.  )  To  insure  that  parking  orbit  points  that  would  produce  a  better  solu¬ 

tion  if  used  are  included  in  the  plane  change  aim  point  arrays. 

Figure  9  shows  a  block  diagram  of  the  plane- change  subroutine  which  follows  screening. 

c.  Ephemeris  of  the  Target  Orbit  The  format  for  expressing  target  orbital  para¬ 
meters  consists  of  a  matrix  of  six  orbital  elements  and  their  first  three  time  derivatives. 
The  parameters  used  are: 

1.  )  Mean  anomaly  (M) 

2.  )  Eccentricity  (e) 

3.  )  Inclination  ( l ) 

4.  )  Right  Ascension  of  Ascending  Node  (ft) 

5.  )  Argument  of  Perigee  (a>) 

6. )  Semi-major  axis  (a) 

7.  )  The  precise  meaning  of  the  elements  is  based  on  their  use  in  the 

Variable  Point  Guidance  and  Targeting  Program. 

Within  the  program,  prediction  of  future  missile  position  for  the  aim  point  ahead  of  present 
position  is  made  on  the  basis  of  a  Keplerian  orbit  starting  with  the  present  position  and 
velocity.  This  orbit  is  known  as  the  osculating  orbit.  If  the  orbital  perturbations,  e.  g. , 
those  due  to  earth  oblateness,  are  appreciable,  the  true  path  will  deviate  from  the  predicted 
path.  If  the  missile  orbit  differs  little  from  the  target  orbit,  as  will  certainly  be  true  in 
the  final  stages  of  the  rendezvous  problem,  the  perturbative  forces  on  the  missile  and 
target  will  be  approximately  equal.  Therefore,  if  the  target  position  is  predicted  on  the 
basis  of  its  current  osculating  elements,  the  relative  error  at  rendezvous  will  be  much 
smaller  than  the  absolute  deviation  in  predicted  position.  In  the  earlier  portions  of  the 
rendezvous  problem,  the  target  and  missile  orbits  will  differ  appreciably  and  small  predic- 
reduced  as  the  two  orbits  are  brought  into  coincidence.  For  these  reasons,  osculating 
elements  are  used  in  the  element  matrix. 

If  the  time  over  which  the  elements  are  to  apply  is  short  compared  to  the  period  of  perigee 
rotation  (which  is  a  minimum  of  eighteen  days),  the  long  period  terms  may  be  expanded  in 
a  power  series  in  time  and  included  with  the  secular  terms.  Otherwise,  the  secular  terms 
alone  will  be  used. 

d.  Earth  Oblateness.  The  derivatives  of  the  elements  due  to  the  principal  harmonic 
term  in  the  expansion  of  the  earth's  field  are  given  in  many  places.  For  reference,  the 
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Figure  9.  Split  P  Block  Diagram 
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article  "Derivation  of  Nodal  Period"  by  Kalil  and  Martikan,  AIAA  Journal,  Volume  1  to 
No.  9,  p  2041  (September  1963)  is  used  here.  Only  secular  terms  are  included. 


where: 


n 
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VGM 

7 


mean  motion  in  osculating  orbit 


Jgm  =  1.  239443589  x  10" 3  sec' 
J  =  |  J2  =  1. 62345  x  10"3 


a  =  semi- major  axis  of  osculating  orbit  in  earth  equatorial  radii 
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The  provision  for  changing  targets  was  implemented  by  setting  aside  a  storage  register  for 
the  time  at  which  a  new  ephemeris  is  to  be  read.  If  no  change  in  target  is  to  be  made  this 
register  is  set  to  a  very  large  number.  When  the  time  of  effectivity  of  a  new  ephemeris 
(as  indicated  by  this  register)  is  exceeded,  the  new  ephemeris  is  read,  inserted  into  the 
target  ephemeris  registers  and  processed  to  develop  a  new  flight  plan.  The  new  ephemeris 
time  register  is  then  reset  to  the  large  number. 

To  relate  this  capability  to  a  realistic  situation  where  a  change  in  target  or  in  the  final  orbit 
of  the  interceptor  is  required  it  would  be  necessary  to  build  into  the  guidance  computer  a 
provision  from  accepting  25  words  of  information  through  a  ground  data  link  or  G-om  some 
source  internal  to  the  vehicle.  The  read-a-new-ephemeris  time  register  could  be  reset  at 
any  time  through  the  mission. 

This  capability  provides  the  following  advantages  to  the  vehicle  once  in  flight: 

1.  )  The  ability  to  receive  better  data  on  a  present  target 

2.  )  The  ability  to  change  targets  completely 
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3.  )  The  ability  to  start  a  rendezvous  mission  from  a  quiescent  in-orbit 

state. 

4.  )  The  ability  to  handle  multiple  targets  (or  to  inject  multiple  payloads) 

starting  rendezvous  n  +  1  after  completing  rendezvous  n 

5.  )  The  ability  to  adapt  to  a  safe  orbit  if  the  rendezvous  mission  is 

aborted 

6.  )  The  anility  to  accept  orbit  change  instructions  from  a  human  pilot. 
The  ephemeris  of  the  target  at  any  time  is  given  by: 


(m-25) 


e.  Ephemeris  Conversion.  The  ephemeris  is  converted  to  position  and  velocity 
vectors  at  any  time  by  the  following  elliptical  orbit  relationships: 
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(Ill- 34) 


IT  =  RIr 

-  cos  ft  sin  (v  +  co)  -  sin  ft  cos  (v  +  &)  cos  i 

-  sin  ft  sin  (v  +  w)  +  cos  ft  cos  ( v  +  o>)cos  i 
cos  (v  +  u>)  sin  i . 


VR  =  VH  tan  y 
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f.  Position  and  Velocity  of  the  Launch  Site.  Before  the  system  is  switched  to  the 
inertial  mode  the  position  and  velocity  of  the  launch  pad  in  inertial  space  is  computed  from 
the  latitude  and  longitude  as  follows: 
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During  the  prelaunch  phase,  the  primary  function  being  performed  by  the  on-board  computer 
is  targeting  and  flight  planning.  The  line  of  nodes  is  defined  90  degree  downrange  from  the 
launch  site  by  the  following  equation: 


^NODE  ~  *HT  X 
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The  first  aim  point  is  then  set  along  the  defined  line  of  nodes  90  degree  downrange  from  the 
launch  site  at  a  specified  altitude.  The  missile  angular  momentum  unit  vector  is  defined  as 
that  instantaneously  required  to  reach  the  first  aim  point: 


3tj 
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The  impulsive  velocity  (scalar)  required  to  reach  the  first  aim  point  at  zero  flight  path  angle 
is  computed,  along  with  the  velocity  of  arrival.  The  time  of  arrival  is  computed,  including 
a  penalty  (input)  for  losses  during  the  first  burn.  The  impulsive  velocity-to-be-gained 
vector  is  computed.  It  is  assumed  that  the  first  and  second  aim  points  lie  on  the  parking 
orbit.  Two  complete  flight  plans  are  computed  (rendezvous  at  both  nodes)  for  each  lift-off 
time,  and  the  ”better"  cf  the  two  is  chosen.  (Following  lift-off,  only  one  flight  plan  is  com¬ 
puted,  corresponding  to  the  ’’better’'  solution  at  lift  off,  except  under  special  circumstances 
such  as  a  major  shift  of  the  line  of  nodes  or  a  new  target  ephemeris. ) 

During  all  non-thrusting  mission  phases  after  lift-off,  the  missile  unit  angular  momentum 
vector  is  computed  from: 
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"he  sine  of  the  out-of-plane  angle  is  computed  from: 


sin  OOP  -  1HT  x  iHM 


(HI-43) 


If  this  is  less  than  the  sine  of  0. 15  degree,  the  out  of  plane  angle  is  set  to  zero  and  the 
target  line  of  apsides  is  substituted  as  the  line  of  nodes: 
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Otherwise,  the  nodal  unit  vector  is  set  from: 
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The  position  and  velocity  vectors  in  earth  centered  coordinates  are  computed  while  on  the 
launch  pad  by  the  geometric  relationship  of  the  known  latitude  and  longitude  of  the  site  on 
the  rotating  earth. 

After  lift-off  the  position  and  velocity  determination  (navigation)  are  computed  by  integra¬ 
tion  of  accelerometer  pulses  derived  from  an  inertial  measuring  unit. 

3.  Guidance  and  Navigation.  Variable  Point  Guidance  and  Targeting  includes  navigation 
guidance  and  control  subroutines  to  implement  the  flight  plans  generated  by  the  targeting 
procedures.  On  the  basis  of  inertial  sensor  inputs,  the  navigation  subroutine  continuously 
maintains  the  position  and  velocity  of  the  vehicle  over  the  oblate  earth  with  respect  to  an 
inertial  coordinate  system.  The  guidance  subroutines  determine  desired  actions  and 
convert  these  to  steering  commands  and  engine  ignition  and  cut-off  commands. 

Guidance  is  conveniently  discussed  relative  to  the  first  burn,  all  later  major  burns  included 
in  the  flight  plan,  and  one  or  two  minor  burns.  First  burn  guidance  begins  at  lift-off,  the 
latter  being  commanded  by  the  launch  director.  After  an  initial  rise  and  a  roll  to  the 
azimuth  determined  by  the  flight  plan,  the  vehicle  begins  an  open- loop  pitch  program  which 
steers  the  vehicle  through  the  sensible  atmosphere.  Closed-loop  guidance  is  initiated 
midway  through  the  first  burn.  It  operates  to  cause  the  vehicle  to  arrive  at  the  next  aim 
point  (whose  location  may  be  fixed  or  variable)  with  the  correct  flight  path  angle.  A 
velocity- to-  be- gained  vector  is  continuously  computed  as  that  which,  if  instantaneously 
achieved,  would  allow  the  vehicle  to  free-fall  to  the  aim  point  on  the  specified  trajectory. 
Steering  drives  two  components  of  this  velocity  vector  to  zero.  Cut-off  is  commanded  when 
the  third  component  passes  through  zero. 

The  essence  of  the  guidance  problem  for  the  on-orbit  burns  is  to  convert  from  an  impulsive 
flight  plan  to  a  finite-burn  flight  plan  accomplishing  the  same  results.  Each  on-orbit 
burn  is  handled  in  the  same  way.  The  start  of  the  burn  is  timed  so  as  to  locate  the  centroid 
of  the  predicted  acceleration- time  profile  at  the  time  of  the  equivalent  impulsive  maneuver. 
That  this  is  an  element  of  a  finite-burn  trajectory  optimally  approximating  an  impulse  is 
described  in  Reference  3.  The  attitude  of  the  vehicle  prior  to  the  start  of  the  burn  directs 
the  nominal  thrust  direction  along  the  required  aV  vector.  Once  the  burn  is  started,  the 
guidance  computes  a  velocity- to- be- gained  as  the  difference  between  the  required  velocity 
(that  which  allows  the  vehicle  to  follow  the  current  flight  plan)  and  the  actual  velocity. 
Steering  and  cut-off  are  based  on  the  velocity-to-be-gained  vector  exactly  as  in  the  first 
burn. 

One  minor  burn,  typically  requiring  no  more  than  a  few  feet  per  second,  is  introduced 
between  nodes  to  take  up  the  accumulated  inaccuracies  from  control  tolerances  or  other 
perturbations.  In  a  vehicle  typical  of  the  near  future,  a  minor  burn  may  be  accomplished 
by  o*her  than  the  main  engine,  perhaps  by  the  attitude  control  system  in  a  translational 
mode.  The  guidance  steering  for  this  mode  would  be  the  same  as  that  described  above. 
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a.  Position  and  Velocity  (Navigation)  After  Start  of  Inertial  Mode, 
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b.  Angular  Momentum  Vector  of  Target. 
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c.  Aim  Point  One.  Prom  launch  site  the  node  will  be  placed 
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The  radius  vector  of  the  first  aim  point  is  then: 
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Angular  momentum  vector  of  guided  vehicle  on  pad  define: 
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d.  Launch  Azimuth.  The  launch  azimuth  is  used  to  define  the  required  roll  maneuver. 
Range  safety  requirements  may  limit  the  allowable  launch  azimuth. 
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where: 
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and 


1N  =  lg  x  1R  (North  vector) 
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e.  Velocity  to  be  Gained.  Having  computed  the  initial  position  vector  Rj  and  the 
radius  vector  of  the  final  point  Rj  the  sine  and  cosine  of  down  range  angle  between  these 
points  is  computed 


40 


c  =  cos  a.  =  Tr.  .  Tr. 


S  =  sinfiy  =|lR.  xTr.| 
sign  of  s  =  sign  of  (T.  x  V  '  1HT 


tany,  =  I1  ‘  R?)  ’  W.  y  j 
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f.  Steering  During  Thrust 
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h.  Time  to  Start  the  Burn  (Two  Stages).  The  calculation  of  time  to  start  the  on-orbit 
burns  cover  the  possibility  that  two  stages  will  be  required,  because  of  the  imminent  fuel- 
depletion  of  one.  The  equations  for  this  are  based  on  a  constant  thrust,  constant  mass-flow 
rate  vehicle.  In  the  present  version,  once  a  burn  is  started,  the  requisite  number  of  stages 
needed  to  complete  the  burn  will  be  properly  called  for. 
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(2)  If  Both  Stages  Complete  the  Burn: 
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i.  Time  to  Cut  Off  Engine.  The  time  to  cut  off  is  based  upon  the  component  of  Vg 
along  the  1/ steering  axis.  The  other  two  components,  Vgt  and  VgTj  are  continually  driven 
to  zero  by  the  pitch  and  yaw  steering  'ommands.  As  Vg|  oiminishes,  the  time  of  zero 
crossing  is  predicted,  taking  into  account  tne  predicted  tail  off  effects  of  the  engine.  A 
discrete  is  issued  at  the  correct  time. 
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SECTION  IV 


RENDEZVOUS  MISSION  SIMULATION  STUDIES 


Task  4  of  Contract  Ar04(695)-946  requires  that  several  typical  rendezvous  missions  be 
simulated  and  their  results  plotted  and  tabulated.  Twelve  simulation  runs  were  made 
against  the  various  targets  specified  by  contract  (Paragraphs  3.4. 1  through  3.4.4)  and 
the  results  obtained  in  these  runs  are  shown  and  discussed  in  this  section. 

The  particular  targets  specified  were  the  following: 

Target  41  (Paragraph  3.4,  1)  "A  satellite  in  a  circular  orbit  at  an  altitude  of  150  nautical 
miles  inclined  35  degrees”. 

Target  42  (Paragraph  3. 4,  2)  ”A  satellite  in  a  circular  orbit  at  1000  nautical  miles  in¬ 
clined  49  degrees". 

Target  43  (Paragraph  3.4.  3)  "A  synchronous  satellite  in  an  orbit  whose  perigee  is  500 
nautical  miles  at  an  inclination  of  65  degrees". 

Target  44  (Paragraph  3. 4. 4)  "A  satellite  in  an  elliptical  polar  orbit  with  a  200  nautical 
mile  altitude  of  perigee  and  a  300  nautical  mile  altitude  of  apogee  and  90  degree 
inclination". 

The  target  definitions  given  above  do  not  completely  specify  the  rendezvous  situation  for 
these  targets  since  no  initial  conditions  are  given.  In  order  to  make  the  runs  required  a 
set  of  initial  conditions  was  selected  for  each  target.  An  attempt  was  made  to  make  this 
selection  on  a  rational  basis  which  will  be  discussed  below  with  the  results  for  each 
target. 

1.  General  Discussions  of  Simulation  and  Cases  Studied.  Prior  to  making  the  simulation 
runs,  extensive  changes  were  made  to  the  Variable  Point  Simulation  Program  to  conform 
to  Task  1  (Paragraph  3. 1)  of  AF04(695)-946.  These  changes  fall  into  the  following  four 
categories: 

a.  Category  1  -  Program  Improvement  and  Cleanup.  The  effort  performed  in  this 
category  was  primarily  in  reducing  the  number~of  subroutines  required  in  the  program. 
This  was  done  by  eliminating  routines  whose  functions  had  been  superseded  and  by  con¬ 
solidating  functions  of  several  routines.  The  number  of  subroutines  (after  all  changes 
discussed  below  had  been  made)  was  changed  from  96  to  78.  During  the  cleanup  process 
about  3000  words  of  required  storage  were  eliminated.  These  savings  were  realized 
even  after  the  addition  of  the  routines  necessary  to  generate  a  data  tape  for  the  SC4020 
plot  output. 
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Also,  in  this  category,  the  simulation  program  was  completely  converted  from  Fortran  II 
to  Fortran  IV.  This  change  required  more  effort  than  simply  examining  the  language  of 
the  program  and  making  appropriate  changes.  The  7090  Fortran  n  system  allows 
division  by  zero  and  square  roots  of  negative  numbers  without  a  machine  stop.  The  7094 
Fortran  IV  system  does  not.  This  led  to  a  review  of  each  routine  during  which  major 
sections  of  routines  were  redesigned.  Fortran  IV  options  such  as  name  list  and  data  were 
also  incorporated.  Results  obtained  indicated  that  the  conversion  was  successful. 

b.  Category  2  -  Incorporation  of  Tangent- Tangent  Transfers  and  Full  Orbit  Phasing. 
This  change  required  a  redesign  of  the  targeting  routines  that  existed  at  the  initiation  of 
the  contract.  The  major  subroutines  resulting  from  this  redesign  were  PLAN,  which  con¬ 
tains  the  targeting  functions  for  the  major  part  of  the  trajectory,  and  GAME  which  con¬ 
tains  the  logic  to  solve  the  time  constraint  for  rendezvous  using  full  orbit  phasing.  The 
changes  made  were  essentially  those  described  in  Reference  9.  No  real  problems  were 
encountered  in  this  phase  of  contract  effort  and  no  failures  of  the  program  in  these  sub¬ 
routines  have  been  encountered.  The  combination  of  PLAN  and  GAME  performs  the 
functions  of  generating  complete  rendezvous  trajectories  from  the  present  location  of  the 
missile  to  an  internally  generated  rendezvous  point  in  the  future.  The  output  of  these  sub¬ 
routines  is  an  array  of  aim  points  on  this  trajectory  through  which  the  missile  must  fly. 
The  data  stored  for  each  of  these  aim  points  include  the  position  vectors,  radii,  velocities 
and  flight  path  angles  entering  and  leaving,  the  plane  change  to  be  executed  and  the  time 

at  which  the  missile  must  pass  through  each.  This  data  fully  defines  the  planned  missile 
trajectory  to  rendezvous. 

c.  Category  3  -  Further  Refinement  in  Accuracy  of  Closure.  The  effort  involved  in 
this  portion  of  Task  1  resulted  in  subroutine  PLANCM  which  is  described  in  Part  2, 
Volume  I  of  this  report.  The  solution  chosen  was  essentially  that  described  in  Paragraph 
4.  1. 2.  2  of  Reference  9. 

Previously  bi-elliptic  transfers  over  a  simulated  oblate  earth  were  used.  The  range  and 
relative  velocity  conditions  observed  at  the  end  of  a  mission  were  marginally  acceptable. 
The  resulting  miss  distance  and  velocity  deviation  was  a  result  of  the  method  used  for 
generating  the  aim  point  parameters.  Aim  points  were  generated  as  vectors  that  were  on 
the  true  target  orbit  accounting  for  oblateness  effects.  In  contrast,  the  missile’s  flight 
path  between  aim  points  was  predicted  on  a  Keplerian  basis  over  a  spherical  earth.  The 
result  was  that  during  the  coast  time  between  aim  points  the  missile  was  acted  on  by  the 
simulated  earth  oblateness  and  consequently  its  real  position  deviated  from  the  desired 
position  (3. 7  nmi  in  one  case)  from  its  predicted  position  as  the  rendezvous  aim  point  was 
approached. 

A  two  part  solution  was  attempted  for  increasing  the  accuracy  and  performance  in  the 
closure  mode  or  terminal  phase  of  the  mission.  The  first  part  of  the  solution  was  very 
effective  while  the  second  appeared  to  cause  difficulty  in  the  terminal  phase  of  the  simu¬ 
lation  runs  against  Target  43.  the  highly  elliptical  synchronous  target,  though  it  was 
effective  against  the  other  three  lower  targets. 

The  first  portion  of  the  solution  consisted  of  changing  the  method  by  which  aim  points  are 
generated.  The  change  from  the  method  described  above  was  to  predict  future  positions 
of  both  vehicles  consistently  (eliminated  using  ephemeris  oblateness  terms  in  computing 
future  target  positions).  The  present  prsition  is  continuously  computed  by  the  corrected 
navigation  equations  but  the  future  positions  are  computed  by  a  more  simple  means. 
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The  updated  orbital  parameters  are  projected  ahead  with  Kepler ian- spherical  earth 
equations  to  determine  the  position  of  the  node  between  the  two  orbit  planes  and  the  con¬ 
ditions  with  which  the  target  will  arrive  at  the  nodal  points.  The  justifying  assumption 
for  this  change  is  that  the  orbits  of  the  target  and  missile  become  very  close  as  both 
approach  the  rendezvous  point  where  each  will  be  equally  affected  by  earth  oblateness.  Thus, 
the  rendezvous  point  is  predicted  by  projecting  the  orbits  of  both  the  target  and  missile  with¬ 
out  considering  oblateness  effects  in  the  projection  process.  The  actual  intercept  aim  point 
may  drift  slightly  from  the  predicted  rendezvous  point,  but  since  the  target  and  missile  are 
affected  in  a  nearlv  identical  manner  by  earth  oblateness,  they  wiil  both  arrive  at  the  drifted 
point  at  the  same  time.  At  that  time  the  error  in  predicting  future  aim  points  will  approach 
zero,  since  after  adaptation  at  the  rendezvous  point  the  missile  and  target  orb<ts  are 
nearly  identical. 

The  second  part  of  the  solution  was  to  use  a  variation  of  the  Phase  A  bi-eiliptic  phasing 
scheme  (see  Reference  9)  as  a  "closure  guidance"  mode.  In  closure  guidance,  which  is 
entered  as  the  missile  approaches  the  last  half  orbit  of  the  phasing  trajectory,  a  new  set 
of  aim  points  is  set  up  for  vernier  adjustments  to  the  trajectory.  These  aim  points  are 
located  at  the  rendezvous  point  (aim  point  no.  5)  and  22. 5,  45,  90  and  180  degrees  uprange 
from  the  rendezvous  point  (aim  point  numbers  4  through  1,  respectively).  This  new  set 
of  aim  points  is  then  used  to  guide  the  missile  through  the  last  stages  of  the  trajectory. 

As  the  missile  approaches  each  of  the  first,  second  or  third  aim  points  the  altitude  of  the 
aim  point  halfway  between  the  missile  and  the  rendezvous  puint  is  adjusted  to  solve  the 
time  constraint  for  rendezvous.  The  missile  then  executes  the  upcoming  burn  aiming  for 
the  adjusted  aim  point. 

The  object  of  this  scheme  is  to  further  improve  the  accuracy  of  arrival  of  the  missile  at 
the  rendezvous  point.  Any  differential  oblateness  errors  that  may  filter  through  the  early 
phases  of  the  mission  are  successively  attenuated  by  the  shorter  downrange  angles  used 
in  trajectory  prediction. 

The  two  methods  fr-r  reducing  the  errors  in  rendezvous  due  to  oblateness  have  worked 
quite  well  in  all  cases  except  the  highly  elliptical  synchronous  target.  As  an  example  for 
comparison  a  table  of  results  of  runs  made  with  the  same  set  of  initial  conditions  against 
Target  42  is  included  as  Table  V  to  show  the  improvement  in  rendezvous  conditions 
achieved.  The  first  two  runs  shown  in  this  table  are  taken  from  data  shown  in  Reference  9. 
The  first  is  made  without  simulation  of  earth  oblateness  while  both  of  the  others  have 
oblateness  effects  included;  the  deviations  in  range  ;  '  -elative  velocity  from  the  spherical 

case  are  shown  in  parenthesis.  In  this  example  the  changes  described  above  have  cut  the 
deviation  in  relative  velocity  by  50  percent  and  practically  eliminated  the  deviation  in  range. 

In  the  case  of  the  highly  elliptical  synchronous  target  orbit,  termination  or  rendezvous  was 
not  achieved  in  making  the  simulation  runs  when  using  oblateness  effects.  The  principle 
symptoms  encountered  were  failures  in  the  closure  mode  targeting  section  when  it  is  first 
used  about  ICO  degrees  from  the  planned  rendezvous  point.  These  failures  were  of  two 
general  types.  The  first  type  was  a  program  stop  due  to  the  computation  of  a  "negative 
semi-major  axis",  an  indication  that  a  hyperbolic  trajectory  is  required  by  the  missile  to 
satisfy  the  rendezvous  time  constraint.  No  capability  for  flight  time  computation  in  the 
hyperbolic  range  has  been  built  into  the  program,  therefore  this  stop  has  been  set  up  for  de¬ 
bugging  purposes.  The  second  type  of  stop  encountered  was  the  "NOGO,  rendezvous  not 
achieved  within  fuel  budget".  This  stop  is  actuated  when  the  predicted  aV  to  finish  the 
mission  has  exceeded  tne  aV  remaining.  Since  all  liftoffs  were  made  with  a  aV  margin  over 
and  above  the  computed  aV  required,  the  error  stop  indicated  that  some  major  change  had 
taken  place.  The  diagnosis  of  these  error  stops  is  discussed  more  fully  in  Paragraph  4. 
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Table  V.  Comparison  of  Runs  Against  Target  42  (Lift-Off  at  8. 33  hrs) 


Flight 

Time 

Total 

AV 

Range  at 

Rendezvous  Burn 
Completion  (nmi) 

Relative  Velocity 
at  Rendezvous 
Completion  (ft/sec) 

Phase  A* 
Spherical  Earth 

4.432 

33,005 

1.01 

5.65 

Phase  A* 

Oblate  Earth 

4.418 

33,927 

4.71  (+4. 70) 

16.98  (+11.33) 

Run  No.  101 

5.416 

33,320 

1.  26  (+0. 24) 

11.62  (+5.36) 

*Run  Under  Contract  AF04(695)-633 


d.  Category  4  -  Use  of  Reference  Velocity  Vectors  for  Reorientation  During  Coast 
Periods,  ThfJchange  in  the  method  of  orienting  the  vehicle  for  upcoming  burns 
is  contained  almost  completely  in  subroutine  GUIDON  described  in  Part  2  of  the  volume. 

Prior  to  this  contract  reorientation  for  on-orbit  burns  was  begun  at  the  time  the  vehicle 
approached  a  specific  time  increment,  225  seconds,  prior  to  the  upcoming  aim  point 
time.  As  an  example  as  the  vehicle  approached  the  n™  aim  point,  whose  computed  time 
of  arrival  was  1200  seconds,  computation  of  the  velocity -to -be -gained  vector  (Vg)  would 
begin  at  975  seconds.  Vg  would  be  based  upon  a  trajectory  from  the  missile  position  at 
that_time,  975  seconds,  to  the  next  aim  point  downrange,  n  +  1.  The  vector  components 
of  Vg  would  then  be  processed  and  attitude  rate  command®  computed  and  issued  to  align 
the  vehicle  for  the  upcoming  burn.  This  scheme  was  effective  in  almost  all  cases.  In 
cases  where  very  long  on-orbit  burns  were  encountered,  the  vehicle  could  not  reorient 
quickly  enough  to  a  proper  attitude  prior  to  the  start  of  the  burn. 

The  approach  used  was  to  begin  the  computation  of  a  VfT  vector  for  the  next  burn 
immediately  after  completion  of  the  previous  burn.  Tlie  parameters  of  missile  arrival 
at  the  next  aim  point  are  computed  and  compared  with  the  stored  parameters  leaving  that 
aim  point  (velocity,  flight  path  angle  and  plane  change)  to  compute  an  MimpulsiveM  Vg 
vector.  This  Vg  is  then  used  for  reorientation  during  coast.  The  calculatic  is  repeated 
every  100  seconds  during  coast  periods  (see  subroutine  CONTROL)  to  update  the  Vg 
vector  whose  magnitude  is  also  used  to  compute  the  time  to  start  the  next  burn.  As  the 
time  to  start  the  burn  is  approached  the  calculation  reverts  to  the  previous  method  of 
computing  Vg  and  vernier  corrections  to  attitude  are  commanded.  Using  the  example 
above,  a  typical  sequence  would  be  to  compute  an  impulsive  Vg  and  time  to  start  a  burn 
at  800  seconds.  These  values  would  be  updated  at  900,  1000  ar»i  1100  seconds  reissuing 
attitude  commands  at  each  computation.  For  a  burn  startingjai  1150  seconds,  switch¬ 
over  would  occur  at  1120  seconds  to  the  previous  method  of  Vg  computation. 

2.  launch  Trade-Off  Charts.  Prior  to  running  the  full  simulation  program  in  most  cases, 
a  ’’launch  trade  -off  s  ranaly  sis  is  usually  made  to  determine  the  parameters  of  the  rendez¬ 
vous  situation.  Program  PHIBOP,  an  impulse  program  which  incorporates  the  essential 
ideas  of  the  flight  planning  part  of  Variable  Point  Guidance  and  Targeting  is  particularly 
suited  to  this  purpose.  The  inputs  to  this  program  are  target  ephemeris,  launch  site 
location  and  time  period  to  be  investigated.  The  program  computes  various  parameters 

of  the  rendezvous  situation  such  as  total  AV,  flight  time,  out-of-plane  angle  and  launch 
azimuth  as  a  function  of  lift-off  time  relative  to  some  epoch  clock  time.  From  these 
parameters  lift-off  times  oi  interest  may  be  chosen  and  the  full  Variable  Point  Guidance 
and  Targeting  simulation  program  used  at  these  lift-off  times  to  obtain  complete  infor¬ 
mation.  Program  PHIBOP  was  modified  prior  to  this  contract  to  be  consistent  with  the 
inputs  for  program  SCPLOT,  the  SC4020  plotting  program  described  in  Section  IV  of  Par1.  2, 
Volume  I.  The  output  data  of  interest  from  PHIBOP  was  plotted  and  used  in  many  cases 
to  select  lift-off  times.  Other  data  useful  as  input  to  the  simulation  program  is  also 
obtainable  from  these  plots,  i.e. ,  from  flight  time  output  a  maximum  simulated  flight 
time  (TFINAL)  may  be  selected, 

3.  Guided  Multi-Stage  Vehicle  and  Environmental  Simulation  Results.  Paragraphs  a. 
through  d.  discuss  the  initial  conditions  used  for  the  various  ephemerides  and  compare  re¬ 
sults  with  previous  results  where  available.  A  tabular  summary  of  results  obtained  in  the 
twelve  contractually  required  demonstration  runs  is  included  as  Table  IV.  All  runs,  except 
those  for  the  synchronous  ellipse,  were  made  with  earth  oblateness  simulated.  These  re¬ 
sults  will  be  discussed  in  the  following  sections. 
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Table  VI.  Variable  Point  Guidance  and  Targeting  Demonstration  Run  Results 


Relative 
Velocity 
Prior  to 
Final  Burn 
(ft/sec) 

1319.4 

452.3 

113.  1 

632.5 

23.8 

486.7 

— 

t-  00 

cm*  *-!  co 
■M*  t-  o> 

H 

189.4 

553.7 

179.3 

WADC 
Output 
Tape  No. 

2630 

2640 

2641 

t-  CO  to 
NC4N 

CO  CD  CD 
CM  CM  CO 

Tr  O)  oi 

co  cm  co 

CO  CD  CD 

CM  CM  CM 

O  CM  CO 

CM  •M*  tT 

CO  CD  CD 
eM  CM  CM 

Flight  Time 
to  100  nmi 
Range  (hrs) 

2.627 

3.402 

2.557 

in  CM 

CM  CM  CO 

H  CM  CO 

in  oo  r- 

r-  co  cc 

O  CM  co 

in  tt  m 
rr  rt  rr 

2.883 

4.976 

4.951 

Relative 
Velocity 
(ft,  sec} 

13.98 
12.76 
12.  54 

0.  84 

8.  €0 
4.90 

j 

CM  CO  CM 

CM  TT  CO 

6  w  rr 

CO  CM  CO 
iroto 

co’  in 

Wt.  In 
Orbit 
(lbs) 

8,  617.7 
21, 000 
27, 000 

20,  600 

8,  000 
21,180.7 

5,  165.6 
5900. 
6500. 

_ 

19,  000 

6,  588.  3 
8,314.5 
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(nmi) 

cm  «o  r- 
f-  to 

CO  CO  CO 

1.40 

1.71 

1.74 

rtOOt^ 

Of-H 
d  O  r-J 

1.  80 
2.00 
1.  15 

Flight 
Time  (hrs) 

2.763 

4.83 

3.833 

5.416 

9.060 

7.998 

46.  73 
46.97 

47.  82 

3.  526 
5.044 
5.055 

Total 

AV  (ft/sec) 

38,  411 
33,011 
30,530 

33,320 

39,264 

32,942 

43,572 

42,463 

41,533 

33,742 

41, 118 
38,785 

Run 

No. 

_ 

r*  CM  CO 
CO  CO  CO 
H 

rt  CM  O) 
o  o  o 

H  H  H 

r-t  CM  CO 

7M  CM  CM 

T-r  r-<  r-t 

^  CM  CO 

T*  r-» 

*— 1  *H  1-4 

Lift- 

Off 

Time 

(hrs) 

in 

OHt» 

co  co  'O’ 

8.25 

3.0 

0.33 

O  O  CO 

in  o  co 

fH  H  o 

CM  C-  CO 

CM*  d  ro 

Angle 
out -of 
Plane 

CO 

HON 

tD  tD  H 

0.2 

18.8 

2.90 

19.73 
15.  00 
8.  00 

5.  70 
23.  83 
14.04 

150  nmi 
Circular 

35  degrees 

1000  nmi 
Circular.. 

49  degrees 

— 

Synchronous* 
Ellipse  at 

65  degrees 

500  nmi 
Perigee 

Polar 

Ellipse 
200-300  nmi 

•on  ia9.rei 

t— 1 

42 

43 

44 

3 

& 

JC 

a 
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O 
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u 
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u 
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o 

s 

3 
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In  all  of  these  runs  a  parking  orbit  altitude  of  100  nmi  was  specified.  Therefore,  in  each 
case,  the  simulated  vehicle  flew  to  a  first  aim  point  at  100  nmi  altitude.  90  degrees  down 
range  from  the  launch  point.  The  following  group  of  plots,  Figure  10  through  16,  represent 
a  time  history  of  the  first  burn  oi  a  simulated  Titan  IIIC  launch  vehicle.  The  particular 
mission  flown  was  a  rendezvous  mission  against  Target  42;  this  run,  number  100,  was 
later  repeated  with  oblateness  to  become  run  number  109,  The  launch  azimuth  involved 
was  48. 8  degrees.  The  parameters  plotted  as  functions  of  time  include: 

1. )  Altitude 

2. )  Velocity 

3. )  Mach  number 

4. )  Dynamic  pressure 

5. )  Body  rate  commands  in  radian  s/second 

6. )  Alpha  -  angle  of  attack  of  pitch 

7. )  Beta  -  angle  of  attack  in  yaw. 

Mach  number,  dynamic  pressure,  alpha  and  beta  are  set  to  zero  after  core  2  cut-off  when 
the  atmospheric  calculations  are  discontinued.  The  plot  of  body  rate  commands  (Figure  14) 
shows  various  phases  of  the  first  burn  very  clearly.  This  plot  is  a  3  ordinate  plot  in  which 
the  star  (*)  represents  pitch  commands  and  the  cross  (X)  roll  commands.  Lift-off  occurs 
at  8  seconds.  The  yaw  command  (represented  by  a  straight  line)  and  roll  commands  arc 
both  zero  throughout  the  mission  except  between  10  and  16  seconds  where  the  roll  program 
is  invoked.  The  commands  for  the  pitch  program  can  be  seen  between  8  and  120  seconds 
and  the  new  pitch  command  sequence  can  be  seen  starting  as  each  new  core  is  initiated. 

One  yaw  command  was  issued  at  486  seconds,  probably  to  provide  a  yaw  vernier  as  the 
transtage  was  started. 

The  dynamic  pressure  and  angle  of  attack  were  found  to  be  well  within  maximum  safe 
values  for  the  Titan  IIIC  specified  in  Paragraph  3. 5. 1  of  Reference  10. 

a.  Target  in  Circular  Orbit  at  150  nmi,  i  =  35  degrees.  The  initial  conditions 
selected  for  this  target  were  based  on  the  hypothetical  case  of  a  satellite  launched  from 
AMR  into  a  35  degree  nclination  orbit.  The  emphemeris  is  written  so  that  at  time  zero 
the  satellite  is  about  20  degrees  downrange  from  its  launch  point.  The  ephemeris  inputs 
computed  for  this  situation  are: 


M  = 

68.3 

M  =  +  5759. 9488 

Q  = 

53.36 

h  =  -  7.0281238 

O)  = 

0.0 

w  =  +  10. 115821 

i  = 

35° 

e  = 

0.0 

a  = 

1.0435548 

The  first  derivative  terms  are  computed  from  Appendix  IV  of  Part  2  of  this  volume.  The 
hypothetical  interceptor  vehicle  is  based  at  AMR  and  the  rendezvous  situation  is  to  be  in- 
vestigate^khrough  a  10  hour  period. 

Program  PHIBOP  was  run  using  this  data  and  used  to  generate  the  six  plots  shown  in 
Figure  17  through  22.  Each  of  these  shows  a  time  history  of  one  item  of  interest  as  a 
function  of  rendezvous  time  throughout  the  10  hour  period.  As  an  example,  a  lift-off  at 
1.  0  hours  would  result  in  a  flight  timr  of  7.6  hours  (Figure  17)  a  total  AV  of  32,200  FPS 
(Figure  18),  an  out-of -plane  angle  of  4.  5  degrees  (Figure  19)  and  a  launch  azimuth  of 


53 


Tt  degrees  (Figure  20).  Figure  21  displays  one  of  the  intermediate  quantities  used  in  the 
calculation,  the  AV  required  for  the  rendezvous  aim  point  burn.  This  quantity  is  equivalent 
to  the  relative  velocity  between  the  missile  and  target  satellite  prior  to  the  rendezvous 
burn  (terminal  maneuver).  Figure  22  displays  another  of  the  intermediate  quantities  used 
showing  the  period  of  the  phasing  orbit  in  hours.  This  value  is  about  1. 5  hours.  The 
comparatively  brief  periods  when  the  plot  flattens  out  below  this  value  indicates  the  times 
when  the  mi5  sile  may  launch  and  stay  completely  below  the  target  orbit  altitude  up  to  the 
rendezvous  point. 


Three  simulation  runs  (131,  132,  133)  were  made  against  this  target.  The  summary  results 
are  shown  in  Table  VI.  The  range  and  relative  velocity  figures  shown  at  the  end  of  the  rendez¬ 
vous  burn  are  about  the  worst  shown  for  any  of  the  four  targets  averaging  about  3. 7  nmi  and 
13  ft/ sec.  All  of  these  runs  were  continued  beyond  the  nominal  rendezvous  point  to  deter¬ 
mine  if  the  range  and  relative  velocity  errors  would  improve.  The  range  and  relative  velocity 
results  of  a  half  orbit  after  the  nominal  rendezvous  burn  are  shown  in  Table  VII  with  reductions 
in  errors  shown  for  comparison.  Range  error  is  not  improved  significantly  but  relative 
velocity  error  is  cut  an  average  of  25  percent  for  the  three  runs  by  the  extra  pass  through 
the  closure  mode  guidance  scheme.  Further  extra  passes  did  not  reduce  the  error  beyond 
this  value  indicating  that  a  steady  state  condition  had  been  reached.  It  is  suspected  that 
input  ephemeris  inconsistencies  may  account  for  some  portion  of  the  range  and  relative 
velocity  errors  as  well  as  affecting  the  closure  guidance  scheme  itself. 

b.  Target  in  Circular  Orbit  at  1000  nmi,  Inclined  49  Degrees.  This  target  had  been 
used  as  a  demonstration  case  in  two  previous  Variable  Point  Guidance  and  Targeting  con¬ 
tracts.  Initial  conditions  used  were  the  same  as  these  used  in  the  previous  Variable  Pemt 
Guidance  and  Targeting  final  reports  (References  9  and  11).  At  time  zero  the  target  is 
placed  at  its  ascending  node  and  the  target  plane  passes  through  the  launch  site  (OOP  =  0) 
at  Cape  Kennedy.  A  PHIBOP  run  was  made  for  this  situation  over  a  10  hour  period  and  the 
results  are  shown  in  Figure  23  through  28.  These  plots  are  read  in  the  same  manner  as 
those  of  Figure  7  through  22. 


M 

=  0.0 

M 

=  4186, 5437  °/day 

Q 

=  -28.  145330 

tl 
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Three  simulation  runs  were  made  to  fulfill  the  demonstration  requirements  of  the  contract. 
The  results  o'  these  runs  are  shown  in  Table  VI.  Ail  runs  were  made  with  earth  oblateness 
simulated.  The  runs  were  made  with  closure  criteria  settings  (range  =  15,000  feet  and 
relative  velocity  =  15  ft/second)  that  caused  the  program  to  stop  at  the  end  of  the  nominal 
rendezvous  burn.  Two  cases  are  mailable  from  Reference  9  for  comparison  with  runs 
made  on  this  target.  The  first  of  these  has  been  previously  shown  in  Paragraph  1  to  illu¬ 
strate  improvement  in  accuracy  of  closure.  It  should  al5*o  be  noted  from  Table  V  that  the 
total  AV  required  for  this  mission  has  been  reduced  by  600  feet  per  second  by  incorporation 
of  the  tangent -tangent  transfer  scheme. 

The  second  comparison  case  available  can  be  made  with  runs  that  did  simulate  oblate  earth 
conditions.  Run  No.  100  (109  without  oblateness)  is  directly  comparable  to  the  run  shown 
in  Table  VI,  Reference  9.  Various  parameters  from  these  runs  are  shown  in  Table  VIII 


Table  VII.  Improvement  in  Rendezvous  Conditions 


(Run  100)  7.66  32,924  0.07  0.85 

Phase  A  1.46  33,972  0.78  5.78 


which  shows  an  improvement  of  over  1000  feet  per  second  in  total  AV  for  'die  mission  at  a 
cost  of  6. 2  hoard  of  flight  time.  These  results  were  better  than  those  expected  at  the 
writing  of  Reference  9  and  the  initiation  of  this  contract. 

Run  100  was  also  used  to  generate  several  samples  of  SC4020  plots  of  the  type  required. 
These  are  included  as  Figure  29  through  32  and  show  time  histories  of  range,  relative 
velocity,  range  rate  and  line  of  sight  from  the  time  range  goes  through  100  nmi  until 
rendezvous.  The  slight  step  functions  in  the  relative  velocity  range  rate  time  histories  at 
26,610  and  27,075  seconds  represent  the  midcourse  correction  burns  at  closure  mode  aim 
points  Number  3  and  4.  The  ramp  between  27, 535  and  27, 565  shows  the  relative  velocity 
being  driven  toward  zero  by  the  rendezvous  burn  of  transtage. 

c .  Target  in  Elliptical  Syn  Chronous  (24  Hour)  Orbit  Inclined  65  Degrees  With  Specified 
Altitude  of  Perigee.  The  initial  conditions  selected  for  this  target  were  based  on  the 
hypothetical  case  of  a  satellite  launched  northerly  from  Tyura-Tam  (45.  6  degrees  N  Lati¬ 
tude,  63.2  degrees  E  Longitude)  into  an  elliptical  orbit  with  its  perigee  about  12  degrees 
downrange  from  the  launch  site.  Using  an  earth  rotation  rate  of  0.  72921155  radians  per 
second,  or  360.  98562  degrees  per  day,  results  in  a  semi- major  axis  of  6.  6107069  earth 
radii.  Using  the  perigee  altitude  of  500  nmi  results  in  an  eccentricity  of  0.82676697  (and 
an  apogee  altitude  of  over  3b,  000  nmi).  Toe  other  inputs  were  selected  so  that  the  target 
v/as  at  perigee  and  the  ascending  node  at  20  degrees  west  of  the  launch  site. 

M  =  0.0  M  =  360.96863  °/dav 

ft  =  -20.0  ft  =  -  0.055031037  °/day 

«  =  0.0  b  =  -  0.0069644605  °/day 

i  =65.0 
e  = .82676697 
a  =6. 6107070 

These  inputs  were  used  to  run  program  PHIBOP  and  to  generate  the  plots  included  as 
Figure  33  through  36. 

Figures  35  and  36  show  an  interesting  situation  that  is  troublesome  in  terms  of  launch 
geometry.  At  the  time  that  the  launch  azimuth  is  within  the  normal  AMR  boundaries  44 
degrees  to  110  degrees  (from  Reference  12)  the  out-of-plane  angle  is  high  (over  22  degrees); 
this  is  caused  simply  by  the  65  degrees  inclination  of  the  target  plane  and  results  in  high 
AV  requirements  for  rendezvous  in  this  lift-off  time  period.  Figure  34  shows  that  the 
shortest  flight  times  are  a’so  within  this  range  of  lift-off  times. 

Three  simulation  runs  were  made  against  this  target  using  a  spherical  earth  simulation. 

The  results  of  these  are  shown  as  runs  121,  122  and  123  of  Table  VI.  An  illustration  of 
the  relative  positions  of  the  rendezvous  points  is  shown  as  Figure  39. 

Several  unsuccessful  attempts  vere  made  to  complete  these  runs  using  the  oblate  earth 
effects  of  the  simulation  program.  In  each  case  program  stops  were  encountered  during 
the  last  180  degrees  of  the  mission.  These  stops  are  built  into  the  simulation  program  to 
avoid  unrealizable  and  undesirable  situations.  The  two  messages  encountered  were 
"Nogo,  Rendezvous  Not  Achieved  With  Fuel  Budget”  and  ”Semi-Major  Axis  Negative”; 
these  are  explained  below. 
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Since  no  provision  has  been  built  into  Variable  Point  Guidance  and  Targeting  for  computing 
fligh  .imes  for  parabolic  and  hyperbolic  conic  segments  an  error  stop  has  been  built  in  to 
terminate  the  program  whenever  such  a  segment  is  encountered.  Up  to  this  time  the 
Variable  Point  Guidance  and  Targeting  Program  had  not  been  required  to  rendezvous  with 
a  target  of  this  tyne,  and  no  acceptable  flight  plans  used  had  involved  hyperbolic  or  para¬ 
bolic  flight  plan  segments.  Since  such  segments  require  high  AV  injections  and  "de¬ 
injections"  it  is  doubtful  that  they  would  be  desirable  in  earth  satellite  rendezvous  missions. 
In  addition  this  particular  stop  has  indicated  that  the  flight  plan  sensitivity  to  deviation  is  too 
high. 

The  fuel  budget  stop  is  included  to  indicate  that  the  predicted  A  V  to  finish  the  mission  has 
exceeded  the  aV  capability  remaining  in  the  interceptor  vehicle  at  that  time.  The  maxi¬ 
mum  fuel  budget  for  the  simulated  rendezvous  mission  is  an  input  to  the  Variable  Point 
Guidance  and  Targeting  Program.  The  aV  required  for  the  rendezvous  mission  is 
reasonably  well  known  from  program  PHIBOP.  The  fuel  budget  is  usually  set  slightly 
above  that  amount.  When  the  fuel  budget  error  stop  is  printed  out,  it  generally  indicates 
that  a  major  change  has  occurred  in  the  flight  plan. 


Investigation  of  the  stops  that  occurred  during  attempts  at  making  oblate  earth  runs  in¬ 
dicated  a  problem  in  the  design  of  the  closure  mode  guidance  technique  for  targets  of  this 
severe  eccentricity.  The  results  of  these  investigations  and  probable  reasons  for  in¬ 
adequate  terminal  mode  targeting  are  discussed  fully  in  Paragraph  4. 

d.  Target  in  Elliptical  Polar  Orbit  With  Altitude  of  Perigee  of  200  Nautical  Miles 
and  Altitude  of  Apogee  of  300  Nautical  Miles.  The  hrtial  conditions  selected  for 
this  target  were  such  that  the  perigee  was  at  the  equator  and  ihe  target  was  at  perigee. 
Figure  39  illustrates  the  initial  geometric  relationship  between  the  target  and  the  launch 
site.  A  six  hour  time  scan  was  set  up  for  PHIBOP  with  the  ascending  node  placed  such 
that  a  series  of  southerly  launches  would  be  obtained  from  the  flight  plan  generator.  The 
input  ephemeris  is  as  follows: 


M  =  0  degrees 

M  =  5522. 5555  °/day 

ft  =  220  degrees 

f2  =  0.0 

w  =  0  degrees 

i  =  3.900739 

i  =  90  degrees 

e  =  0.01353572 

a  =  1.0725914  earth  radii 

The  plots  resulting  from  this  PHIBOP  run  are  shown  in  Figure  41  through  45.  Three  lift¬ 
off  times  were  selected  for  simulation  and  the  results  obtained  are  listed  in  Table  VI  under 
run  numbers  111,  112  and  113. 

The  launch  site  used  in  these  runs  was  PMR  (34.  8  degrees  N.  Latitude).  The  only  change 
required  to  the  simulation  program  for  this  launch  site  was  to  change  the  launch  azimuth 
computation.  The  computation  was  changed  to  quadrant  azimu;h  from  0  degrees  to  J6Q 
degrees  rather  than  from  -180  degrees  to  +180  degrees.  The  change  was  made  to  insure 
that  the  launch  vehicle  would  roll  in  the  shortest  direction  during  the  vertical  rise  "roll  to 
azimuth"  maneuver.  This  change  will  not  affect  simulation  of  AMR  launches. 

This  target  was  also  used  to  generate  the  terminal  phase  plots  shown  in  the  folk  ;ng 
section. 
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4.  Terminal  Conditions.  The  following  grouping  of  SC'4020  plots,  Figures  41  through  59 
are  included  to  show  compliance  to  Paragraphs  3. 1.4  and  3.  2  of  contract  AF04(695)-946. 
The  plots  were  made  using  program  SC  PLOT  with  the  output  tape  of  run  111. 


The  plots  cover  the  time  span  in  the  flight  during  which  relative  range  between  the  vehicle 
and  target  was  less  than  100  nmi.  Table  IX  is  a  tabular  representation  of  the  major  events 
that  occurred  within  that  time  span  for  run  111. 

In  run  111  Figures  33  through  37  are  plots  of  the  time  histories  of  the  various  position  and 
velocity  vectors  of  the  target  and  chase  vehicle  as  well  as  the  relative  position  (range)  and 
relative  velocity  between  them.  Since  the  simulation  program  operates  with  different 
integration  time  intervals  in  coast  periods  (10  seconds)  and  burn  periods  (2  seconds)  the 
times  when  the  burns  are  occurring  are  readily  recognizable  from  these  plots  by  the  point 
grouping. 

Figure  38  shows  the  body  rate  reorientation  commands  being  issued  to  the  chase  vehicle. 
The  comparatively  large  ’’humps”  in  pitch  and  yaw  occurring  after  each  burn  illustrate 
the  missile's  reorientation  to  the  reference  velocity  vector  for  the  next  burn.  This  effect 
is  also  apparent  in  a  plot  of  range  vector  in  body  coordinates  included  as  Figure  54. 

Range  rate  time  history  is  shown  in  Figure  55. 

The  next  two  plots  56  and  57  show  the  relative  angular  position  of  the  missile  in  target 
centered  coordinates.  The  center  axis  for  this  coordinate  system  is  the  target  forward 
horizontal  direction.  As  an  example  at  11,400  seconds  the  missile  is  55  degrees  to  the 
■dght  and  45  degrees  below  the  target  forward  direction. 

^  re  58  is  a  multiple  plot  of  line  of  sight  rate  and  its  horizontal  and  vertical  components, 
i*  magnitude  and  vertical  component  appear  equivalent  while  the  horizontal  component 
is  very  near  zero.  Note  that  each  burn  drives  line  of  sight  rate  down  significantly  and 
that  the  final  burn  drives  it  from  its  maximum  value  of  about  2  degrees  per  second  to  zero. 

Figure  59  is  an  example  of  one  of  the  many  cross  plot  combinations  available  using 
SCPLOT  with  the  Variable  Point  Guidance  and  Targeting  simulation  program  output.  A 
listing  of  the  68  quantities  available  for  plotting  is  shown  in  Table  III  of  Appendix  IV, 

Part  2. 

5.  Oblateness  Effects.  During  the  two  previous  contracts  in  which  Variable  Point 
Guidance  and  Targeting  was  developed,  little  emphasis  was  placed  on  making  simulation 
runs  with  earth  oblateness  effects.  The  few  runs  that  were  made  demonstrated  that  the 
simulation  worked  correctly  and  showed  range  and  relative  velocity  errors  at  rendezvous 
that  were  greater  than  those  encountered  in  spherical  earth  runs  but  were  considered  as 
ma:  ginally  acceptable.  A  more  thorough  look  at  the  effects  of  oblateness  on  Variable  Point 
Guidance  and  Targeting  was  initiated  to  improve  the  overall  accuracy.  In  addition  the  re¬ 
quirement  for  full  orbit  phasing  introduced  longer  nominal  coast  periods  between  burns  and 
this  would  tend  to  increase  the  mi^s  distance.  All  runs  shown  in  Table  VI  were  attempted 
with  earth  oblateness  in  effect.  Not  all  of  these  runs  were  successful  in  the  terminal  phase. 
Because  of  the  effect  of  earth  oblateness  in  combinations  with  the  high  eccentricity  and 
period  of  target  43  the  oblate  earth  runs  against  this  target  did  not  terminate  due  to  the  high 
velocity  requirement  generated  by  the  closure  mode  guidance  and  were  rerun  using  a 
spherical  earth  model.  The  error  stops  encountered  were  discussed  in  more  detail  in 
Paragraph  1  and  2.c  and  the  conditions  under  which  they  occurred  are  described  below. 


*•  •  ••  ■■ 
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A  typical  situation  for  the  type  of  failues  encountered  was  the  first  entry  into  closure  mode 
guidance  PLANCM.  This  was  evi  ^nced  by  approximately  five  unsuccessful  attempts  to 
terminate  ’’oblate  earth”  rendezvous  missions  against  the  synchronous  elliptical  target  orbit. 

Since  Target  43  has  such  a  large  period,  phasing  orbits  within  it  will  also  have  large 
periods  leading  to  long  coast  periods  between  aim  points.  During  these  long  coast  periods 
the  orbits  tend  to  drift  slightly  v/ith  respect  to  each  other  in  space.  The  high  eccentricity 
of  Target  43  compounds  the  effect  of  drift  by  introducing  a  secondary  effect.  As  the  orbit 
planes  move  relative  to  each  other,  the  nodal  intersection  line  moves  along  the  orbits  in 
a  toggle-type  motion.  Because  the  nodal  intersection  is  at  a  point  on  the  target  orbit 
different  than  that  originally  designated  for  rendezvous  the  radius,  velocity  and  flight  path 
angle  are  also  changing.  PLANCM  *  erupts  to  fit  a  new  trajectory  through  the  aim  points 
which  it  has  set  up.  This  trajectory,  which  solves  for  rendezvous  under  the  new  nodal  con¬ 
ditions,  includes  an  aim  point  at  an  adjust  a  altitude  to  form  a  two  segment  flight  plan 
rather  than  the  single  segment  flight  plan.  The  use  of  a  large  constrained  flight  path  angle 
at  rendezvous  tends  to  result  in  a  two  segment  flight  plan  that  is  different  enough  in  shape 
from  the  original  phasing  orbit  to  require  excessive  velocity  change.  This  difference  in 
shape  always  results  in  a  total  A  V  larger  than  that  originally  planned;  the  added  aV  was 
large  enough  to  cause  activation  of  the  fuel  budget  error  stop.  When  the  shape  change  was 
more  pronounced  it  indicated  that  an  hyperbolic  type  trajectory  was  required  activating  the 
semi-major  axis  error  stop. 


Table  DC.  Final  Phase  of  R’in  111 


Time  (sec) 

Range  (nmi) 

Event 

10,379.5 

100 

Range  less  than  100  nmi 

11,249.5 

34.4 

Start  transtage  (16,000  lb  thrust)  for  midcourse  correction 
at  closure  mode  aim  point  2  (90  degrees  to  go) 

11,259.2 

33.8 

Cut-off  (aV  was  244.  1  fps) 

11,966.3 

17. 1 

Start  transtage  for  midcourse  correction  at  closure  mode 
aim  point  3  (45  degrees  to  go) 

11,968.2 

17.  1 

Cut-off  (aV  was  45.4  fps) 

1?  ^26.2 

12.  1 

Start  transtage  for  midcourse  correction  at  closure  mode  aim 
point  4  (22.5  degrees  to  go) 

12,329.5 

12.1 

Cut-off  (aV  was  102.3  fps) 

12,685.5 

1.90 

Start  transtage  for  rendezvous  burn  (thrust  =  16,000  lbs, 
total  weight  -  19, 941) 

12,692.7 

l.bu 
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The  situations  encountered  show  a  high  sensitivity  in  the  closure  mode  to  condition  changes 
produced  by  oblateness  in  rendezvous  runs  against  Target  43.  This  high  sensitivity  is  con¬ 
sidered  a  fault  in  the  engineering  design  of  the  closure  mode  PLANCM  .  A  recommend¬ 
ation  for  the  improvement  of  PLANCM  is  included  in  the  following  section. 

An  investigation  was  made  to  determine  whether  or  not  the  methods  of  computing  the  ab¬ 
solute  locations  of  the  target  satellite  and  chase  vehicle  were  compatible.  In  the  simulation 
program,  the  chase  vehicle  position  and  velocity  are  determined  by  integration  of  the 
equations  of  motion.  The  target  position  and  velocity  are  found  by  a  Taylor  series  ex¬ 
pansion  of  the  input  target  ephemeris  in  UPDATE  and  mathematical  conversion  of  that 
ephemeris  in  GEOM  to  the  required  vectors  and  orbital  parameters.  The  two  methods 
were  compared  and  a  more  exact  and  consistent  procedure  for  generating  oblate  earth 
first  derivative  terms  was  derived.  The  refined  ephemeris  terms  resulting  from  using 
this  procedure  with  Target  43  were  used  as  input  to  the  simulation.  The  results  of  these 
runs  on  the  synchronous  ellipse  gave  further  indication  that  the  closure  mode  requires 
further  improvement  or  at  least  a  simple  back  up  mode  in  order  to  assure  that  any  general 
rendezvous  mission  will  be  successfully  and  reliably  completed. 

6.  Conclusions  Derived  from  Simulation  Results.  The  results  included  in  Section  4  show 
the  present  version  of  Variable  Point  Guidance  and  Targeting  to  be  more  effective,  efficient 
and  flexible  than  the  version  that  existed  at  the  completion  of  the  previous  contract, 
AF04(695)-633.  The  improvements  made  to  Variable  Point  Guidance  and  Targeting  were  in 
the  areas  of  fuel  efficiency  (evidenced  by  lower  aV  requirements),  lower  computer  storage 
requirements  and  better  closure  conditions. 

The  twelve  simulation  runs  required  by  contract  AF04(695)-946,  Paragraph  3.4,  have  been 
successfully  completed  and  demonstrate  the  effectiveness  of  Variable  Point  Guidance 
against  a  variety  of  targets.  Limitations  to  the  present  design  of  the  terminal  mode  (last 
half  orbit)  guidance  laws  were  discovered  in  the  case  of  extremely  eccentric  targets. 

Seveial  solutions  to  the  limitations  encountered  are  possible.  A  simple  solution  is  to  con¬ 
tinue  with  the  prior  working  plan  until  sensor  lock-on  occurs.  From  here  better  data  is 
available  for  the  closure  mode.  The  use  of  the  two  point  explicit  guidance  solution  for  the 
closure  burns  wrhere  trajectory  shaping  by  flight  path  angle  iteration  would  be  employed  tv 
solve  the  time  constraint  for  rendezvous  is  another  possibly  better  way  to  terminate  the 
mission.  Revising  the  equations  for  such  a  scheme  to  include  valid  hyperbolic  trajectory 
solutions  would  extend  the  capability  of  Variable  Point  Guidance  and  Targeting  to  more 
effectively  rendezvous  with  targets  of  extreme  eccentricity.  Another  alternative  would  be 
to  return  to  the  primary  targeting  mode  and  updating  the  flight  plan  by  slipping  the  rendez¬ 
vous  point  further  downrange  until  sensor  lock-on  occuis. 
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AL  T  I  TUDE  -  NMILES 


Figure  10.  Altitude  During  First  Burn  as  a  Function  of  Time 
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MISSILE  VELOCITY  -  RUN  100 


TIME  -  SECONDS 


Figure  11.  Velocity  During  First  Burn  as  a  Function  of  Time 
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MACH  NUMBER  -  RuN  100 
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Figure  12.  Mach  Number  as  a  Function  of  Time 
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DYNAMIC  PRESSURE 


Figure  13.  Dynamic  Pressure  as  a  Function  of  Time 

64 


66 


BETA  -  DEGREES 


BETA  -  RUN  100 


Figure  16.  Angle  of  Attack  in  Yaw  as  a  Function  of  Time 
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TIME  VS  FLIGHT  T  I  ME  -  RUN  41  150  NMI  CIRCULAR  TARGET  ORBIT 
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Figure  17,  Time  of  Flight  as  a  Function  of  Lift-Off  Time 
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Figure  20.  Launch  Azimuth  as  a  Function  of  Lift-Off  Time 
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DELTA  V  -  RUN  42  -  1000  NMI  CIRCULAR  TARGET  ORBIT 

**oos 


Figure  26.  Total  Velocity  Change  as  a  Function  of  Lift-Off  Time 
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Figure  27.  Launch  Azimuth  as  a  Function  of  Lift-Off  Time 
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RENDEZVOUS  BURN  -  RUN  42 
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Figure  28.  Velocity  Change  in  the 


i  -  1000  NMI  CIRCULAR  TARGET  ORBIT 


mdezvous  Burn  as  a  Function  of  Lift-Off  Time 
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Figure  30.  Relative  Velocity  as  a  Function  of  Time 
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Figure  31.  Range  Rate  as  a  Function  of  Time 
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Figure  22.  Line  of  Sight  Ratio  as  a  Function  of  Time 
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Figure  34.  Flight  Time  as  a  Function  of  Lift-Off  Time 
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Figure  35.  Angle  Out  of  Plane  as  a  Function  of  Lift-Off  Time 
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Azimuth  as  a  Function  of  Lift-Off  Time 


Figure  38.  °hasing  Orbit  Period  as  a  Function  of  Lift-Off  Time 
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Figure  39.  Geometry  of  the  Synchronous  Elliptical  Orbit 
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Figure  40.  Geometry  of  the  Elliptical  Polar  Satellite  Orbit 
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Figure  42,  Total  Velocity  Change  as  a  Function  of  Lift-Off  Time 
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Figure  43.  Out  cf  Plane  as  a  Function  of  Lift-Off  Time 


Figure  44.  Launch  Azimuth  as  a  Function  of  Lift-Off  Time 
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Figure  45.  PKasing  Orbit  Period  as  a  Function  of  Lift-Off  Time 
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Figure  47.  Target  Position  Vector  as  a  Function  of  Time 
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Figure  48.  Missile  Position  Vector  as  a  Function  of  Time 
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Figure  49.  Range  (ECI)  as  a  Function  of  Time 
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Figure  50.  Target  Velocity  Vector  (ECI)  as  a  Function  of  Time 
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Figure  51.  Missile  Velocity  Vector  (ECI)  as  a  Function  of  Time 
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Figure  52.  Relative  Velocity  (ECI)  as  a  Function  of  Time 
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Figure  53.  Body  Rate  Commands  as  a  Function  of  Time 
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Figure  54.  Range  as  a  Function  of  Time 
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Figure  56.  Line  of  Sight  Heading  as  a  Function  of  Time 
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Figure  57.  Line  of  Sight  Elevation  as  a  Function  of  Time 
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Figure  58.  Line  of  Sight  Rate  and  Two  Components  as  a  Function  of  Time 
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Figure  59.  Line  of  Sight  Hate  as  a  Function  of  R  m  ge 
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Figure  60.  Specific  Force  as  a  Function  of  Range 


SECTION  V 


VARIABLE  POINT  APPLICATIONS 


1.  Introduction.  The  analysis  and  simulation  results  on  Variable  Point  Guidance  and  Tar¬ 
geting  Program  have  shown  that  a  guidance  and  targeting  method  has  been  developed  which 
covers  a  wide  variety  of  space  missions.  The  software  characteristics  in  the  present  state 
of  development  can  provide  both  efficient  flight  plans  and  guidance  commands  within  one 
second  from  insertion  of  mission  objective  to  potential  lift-off  time  from  lift-off  through 
ascent,  mid  course  to  terminal  phase  of  the  following  typical  missions: 

Satellite  Rendezvous  (Manned  or  Unmanned) 

1. )  Space  station  resupply 

2. )  Emergency  rescue 

Single  or  Multiple  Objectives  for  a  Single  Launch 

1. )  Entry  into  a  circular  or  an  elliptical  parking  orbit 

2. )  Orbit  transfer 

3. )  Orbit  plane  change 

Placement  or  Replacement  of  Payloads  into  Precise  Orbits 

1. )  Communications  satellites 

2. )  Synchronous,  equatorial,  polar  and  orbits  of  any  inclination  or 

eccentricity 

2.  Satellite  Rendezvous  (Manned  or  Unmanned).  The  simulation  results  of  Section  IV 
indicated  the  general  capability  for  the  Variable  Point  Guidance  and  Targeting  to  provide 
the  launch  operations  trade-offs  as  well  as  efficient  guidance  for  rendezvous  with  targets 
in  extremely  different  orbits. 

The  targets  used  for  demonstration  of  the  capability  of  the  system  were  intentionally 
chosen  to  be  widely  different  and  yet  typical  of  some  projected  missions.  The  list  of 
targets  is  given  in  Section  rv,  Table  V. 

The  boost  vehicle  used  in  the  simulation  was  representative  of  boosters  currently  avail¬ 
able  using  both  AMR  and  PMR  as  a  launch  site. 

a.  Space  Station  Resupply.  The  space  stations  contemplated  for  the  near  future  are 
in  general  of  low  altitude  (i.  e. ,  MOL,  150  nautical  miles).  Future  space  stations  may  be 
higher. 
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Variable  Point  Guidance  and  Targeting,  in  its  present  form,  is  ideally  suited  to  the  task 
of  space  station  re -supply  since  both  efficient  guidance  and  control  signals  and  flight  plans 
are  continuously  available.  The  computer  output  can  provide  rendezvous  for  any  target 
situation.  The  same  methods  are  recommended  to  be  used  for  any  vehicle,  manned  or 
unmanned.  The  results  have  shown  that  if  the  station  is  in  circular  or  elliptical  orbits 
the  method  is  equally  applicable.  The  launch  windows  available  are  only  dependent  upon 
vehicle  capability.  Variable  Point  Guidance  and  Targeting  give  maximum  launch 
windows  efficient  trajectories,  early  rendezvous,  and  flexibility  for  responding  quickly  to 
payload  changes  and  to  non-standard  target  orbit  characteristics. 

The  method  of  solution  provides  the  flexibility  to  handle  alternate  missions  so  that  as  the 
number  of  possible  stations  in  orbit  or  the  number  and  frequency  of  missions  increases, 
the  need  for  more  flexibility  increases.  In  addition,  if  rendezvous  is  required  when  con¬ 
straints  such  as  sun  line  avoidance,  launch  azimuth  restrictions  or  unavoidable  delays  are 
encountered  the  system  flexibility  aids  in  the  probability  of  accomplishing  the  mission  on 
time  since  a  maximum  launch  window  exists  prior  to  the  insertion  of  other  constraints  or 
the  delays  peculiar  to  a  particular  mission  actually  occurs. 

b.  Emergency  Rescue.  A  versatile  method  such  as  the  Variable  Point  Guidance  and 
Targeting  system  becomes  even  more  necessary  if  minimum  reaction  time  is  important. 

The  emergency  rescue  or  resupply  problem  has  been  considered  by  many  workers  as  an 
eventual  necessity.  Variable  Point  Guidance  and  Targeting  can  handle  any  target  orbit 
in  any  orientation.  In  one  computation  cycle,  it  produces  a  flight  plan  which  can  be 
immediately  used  for  guidance  of  the  rescue  vehicle.  This  method  or  its  equivalent  is 
required  for  the  rapid  response  to  the  unplanned  situation  that  the  term  "rescue"  implies. 

It  has  been  shown  that  if  there  is  doubt  about  the  ephemeris  of  the  vehicle  to  be  inter¬ 
cepted  and  rendezvous  is  to  be  accomplished,  a  vehicle  with  Variable  Point  Guidance  and 
Targeting  can  be  launched,  and  an  updated  ephemeris  sent  up  later  when  available. 

An  IR  and  D  study  was  made  by  RCA  (Reference  13)  on  the  "Cost  Effectiveness  Study  of 
MOL  Resupply  and  Rescue".  An  earlier  version  of  the  Variable  Point  Flight  Plan 
Generator  portion  of  the  program  was  used  to  aid  in  the  generation  of  tradeoff  curves  as 
required  for  this  effort.  Three  boosters  were  compared  from  the  points  of  view  of 
economy,  performance,  and  reaction  time.  The  results  indicated  that  the  Titan  IIIC  was 
best  suited  among  the  three  vehicles  considered  for  the  rescue  mission.  The  other  vehicles 
considered  were  Atlas- Agena  D  and  Saturn  IB. 

Present  operational  guidance  and  targeting  methods  such  as  Gemini  are  not  suitable  for 
the  unplanned  emergency  rescue  because  the  planning  for  rendezvous  into  new  orbits  re¬ 
quires  several  weeks  of  effort. 

The  time  spent  in  the  planning  phase  can  be  costly  ^nd  for  unplanned  situations  the  required 
reaction  time  would  not  generally  meet  the  needs  oi  the  resupply  or  rescue  situation  in¬ 
volved.  The  RCA  Variable  Point  Guidance  and  Targeting  system  not  only  can  meet  fast 
reaction  time  requirements,  when  it  is  needed,  but  it  accomplishes  the  cost  saving  function 
of  planning  for  any  rendezvous  mission  even  if  the  time  to  plan  were  available.  As 
pointed  out  previously,  the  targeting,  flight  plan  generator,  and  guidance  are  contained  in 
a  single  flight-borne  computer  package. 

3.  Single  or  Multiple  Objectives  for  a  Single  Launch.  The  Variable  Point  Guidance  and 
Targeting  method  can  be  used  to  place  any  number  of  vehicles  in  any  number  of  desired 
orbits.  This  means  that  one  or  more  payloads  may  be  mounted  on  a  single  booster,  and 
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each  of  these  proceeds  to  its  required  orbital  destination.  Alternately,  a  single  vehicle 
may  be  placed  sequentially  into  a  series  of  desired  orbits.  It  is  necessary  that  specified 
end  condition  and  navigation  information  on  the  vehicle’s  own  position,  velocity  and  time 
are  known,  or  can  be  obtained. 

a.  Entry  Into  a  Circular  or  Elliptical  Parking  Orbit.  The  entry  into  a  specific  orbit 
is  equivalent  to  rendezvous  with  a  pseudo  target  in  the  desired  orbit.  In  the  present  pro¬ 
gram  the  phasing  is  automatically  bypassed  if  there  is  not  a  target  in  the  required  orbit. 

For  synchronous  satellite  orbits  of  any  eccentricity  and  various  inclinations,  nodal 
crossing,  and  argument  of  perigee,  the  desired  target  ephemeris  is  inserted  as  in  any 
rendezvous  problem  which  depends  upon  real  target  and  precise  knowledge  of  its  time 
position  velocity  relationship.  The  recommended  guidance  and  targeting  for  these  missions 
is  the  standard  rendezvous  method. 

4.  The  Flight  Plan  Generator  as  a  Launch  Operations  Aid.  As  pointed  out  in  Section  Til, 
the  Variable  Point  Guidance  and  Targeting  software  performs  the  following  functions: 

1. )  Target  ephemeris  conversion 

2. )  Flight  plan  generation  (targeting) 

3. )  Guidance  and  control  commands. 

Functions  1.)  and  2.)  can  be  performed  by  a  portion  of  the  total  Variable  Point  Guidance 
and  Targeting  system.  This  portion,  called  the  flight  plan  generator,  can  provide  launch 
operations  with  trade-off  data. 

There  are  at  present  chree  small  programs  based  on  the  Variable  Point  Guidance  and 
Targeting  subroutines  which  were  developed  to  aid  in  checking  out  new  ideas,  and  also 
to  define  flight  plans  for  various  contemplated  space  missions.  These  programs  can  be 
run  on  both  the  RCA  355  and  the  IBM  7090  and  7094  computers.  In  addition,  the  SC4020 
plot  routine  has  been  tied  in  so  that  tradeoff  charts  are  readily  produced. 

These  programs  have  been  proved  to  be  useful  in  their  own  right.  Specifically,  they  provide 
the  launch  trade-offs  data  that  is  so  useful  to  a  launch  commander.  The  launch  trade-off 
data  is  presented  as  a  function  of  time  of  lift-off.  From  the  above  mentioned  programs, 
from  one  to  four  possible  choices  of  rendezvous  trajectories  can  be  presented,  involving 
various  fuel  usages  and  mission  times.  This  data  can  be  matched  to  other  constraints  such 
as  range  safety,  readout  time  and  location,  fuel  budget,  etc.  The  data  available  as  outputs 
of  these  programs  are  as  follows: 

1.  )  Total  propellant  required 

2.  )  Time  of  flight 

3. )  Time  of  mission  completion 

4.  )  Launch  azimuth 

5.  )  Total  plane  change  angle  required 

6.  )  Location  and  time  of  each  burn 

7.  )  Number  of  engine  restarts 

8.  )  Change  in  velocity  and  amount  of  plane  change  at  each  burn 

9.  )  Total  launch  window. 
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The  data  is  useful  to  answer  such  questions  as  the  following: 

1. )  Earliest  possible  lift-off  time 

2. )  Latest  time  of  lift-off  for  successful  mission  completion 

3. )  The  time  of  lift-off  for  minimum  fuel 

4. )  Time  of  lift-off  for  a  particular  launch  azimuth 

5. )  Time  of  lift-off  for  minimum  time  to  mission  completion. 

In  an  operational  situation,  this  data  is  computed  in  real  time  by  subroutines  of  the  Variable 
Point  Guidance  and  Targeting  flight  plan  com  stations  in  the  vehicle  guidance  computer. 

It  may  be  desirable  to  have  a  separate  special  or  general  purpose  ground  computer  to  com¬ 
pute  a  properly  integrated  combination  of  the  above  subroutines  for  the  purpose  of  pro¬ 
viding  the  launch  commanders  with  the  trade-off  data  for  a  12  or  24  hour  look  into  the 
future.  This  will  aid  him  in  selecting  the  best  time  of  launch,  scheduling  operations,  ana 
making  rapid  changes  to  plans  in  case  of  unexpected  delays.  At  present  the  above  in¬ 
formation  is  printed  in  tabular  form  or  as  SC4020  plots.  Automatic  display  equipment 
could  be  used  to  present  the  above  data  as  fast  as  it  is  computed. 
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SEC  TION  VI 


RECOMMENDATIONS  FOR  FUTURE  WORK 


1.  Introduction.  The  development  of  the  Variable  Point  Guidance  and  Targeting  Program 
has  progressed  by  a  series  of  evolutionary  changes.  The  program  in  its  present  state 
could  be  used  for  guidance  and  targeting  from  launch  through  midcourse  to  the  terminal 
phase  in  a  typical  rendezvous  mission.  There  are,  however,  several  areas  where  further 
operational  improvements  should  be  made  before  the  program  is  adapted  to  a  specific 
rendezvous  application.  See  Table  X  for  recommended  follow-on  schedule. 

There  arc  also  many  non- rendezvous  missions  in  which  rapid  flight  plans  would  be  useful. 
The  requirement  to  overfly  specific  points  on  the  surface  of  the  earth  is  needed  in  most 
missions  for  the  purpose  of  maintaining  communication  through  telemetry  stations. 

Range  safety  must  be  considered  in  new  missions. 

The  implementation  of  the  Variable  Point  Guidance  and  Targeting  technique  to  specific 
applications  requires  consideration  of  all  the  detailed  hardware  and  system  integration 
constraints. 

General  improvements  in  the  equations  and  logic  which  would  increase  the  flexibility  and 
efficiency  of  the  system  and  extend  the  applicability  of  the  program  are  discussed  in 
Paragraph  2. 

a.  Discussion.  An  optimized  single  burn  ascent  into  a  low  altitude  parking  orbit 
should  be  included  in  the  program.  Consideration  of  yaw  turns  should  be  included  to  in¬ 
crease  efficiency  in  azimuth  limited  launches.  Efficiency  can  also  be  increased  by  com- 
tining  yaw  turns  with  optimal  placement  of  the  line  of  nodes. 

There  are  available  several  simple  and  close  to  optimal  solutions  to  the  single  burn  to  orbit 
requirement. 

These  methods  should  be  examined  and  compared  to  the  more  sophisticated  methods.  One 
method  should  be  integrated  into  the  simulation  program. 

b.  Automatic  Range  Safety  Implementation  and  Simulation.  The  purpose  of  this  task 
is  to  implement  Range  and  Mission  Safety  requirements  into  the  Variable  Point  Targeting 
System  so  that  it  automatically  seeks  a  flight  plan  which  avoids  constraints.  The  items 
to  be  included  in  safety  considerations  are  constraints  on  launch  azimuth  and  on  impact 
point  of  various  stages  of  the  boost  vehicle. 

Ea  ly  integration  of  range  and  mission  safety  into  the  program  will  allow  continuous  study 
of  methods  to  maximize  performance  and  launch  opportunities  in  the  constrained  situation. 
Means  will  be  examined  for  utilizing  off-set  points  and  dog- leg  maneuvers  as  reasonably 
efficient  means  of  working  around  the  constraints. 
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Table  X.  Recommended  Follow-On  Schedule 


c.  Flight  Plans  for  Overflight  of  One  or  More  Specific  vSites.  All  space  missions  up 
to  the  present  time  have  been  planned  to  overfly  telemetry  stations  and  for  monitoring 
system  performance  during  periods  of  thrust,  and  during  other  operational  sequences. 

Flight  planning  for  overflight  of  points  on  the  surface  of  the  earth  should  consider  range 
safety  limitations. 

d.  Terminal  Guidance  Studies.  The  present  variable  point  scheme  for  accent  oper¬ 
ates  on  a  target  ephemeris  input  as  its  basic  knowledge  of  target  data.  The  scheme  can  be 
carried  into  the  terminal  phase,  achieving  rendezvous  to  the  accuracy  of  the  target  ephe¬ 
meris  itself.  In  general,  this  is  not  sufficiently  accurate  for  typical  missions.  Rendezvous 
accuracy  can  be  improved  with  essentially  no  change  to  the  variable  point  formulation 

if  a  means  existed  for  obtaining  a  more  accurate  ephemeris  than  that  which  was  used  to 
control  ascent.  The  radar  or  other  sensor  measurements  on  the  target  in  the  terminal 
phase  provide  the  necessary  information  for  improving  the  ephemeris. 

It  is  proposed  to  develop  equations  and  incorporate  them  into  the  variable  point  formu¬ 
lation  which  will  convert  radar  measurements  on  the  target  to  an  improved  target  ephe¬ 
meris.  An  auxiliary  problem  is  that  of  smoothing,  that  is,  using  redundant  data  acquired 
over  a  period  of  time  to  reduce  the  statistical  error  of  the  prediction. 

e.  Minimum  Reaction  Time  to  Mission  Completion.  The  intent  of  this  task  is  to 
extend  the  VPG  program  to  include  the  minimum  reaction  time  rendezvous  trajectory 
compatible  with  the  full  utilization  of  a  prescribed  propellant  budget,  extending  the  capa¬ 
bility  still  further  tov.ard  the  emergency  rescue  mission. 

There  are  generally  times  in  a  typical  launch  window  when  a  short  range  direct  ascent 
(two- burn)  rendezvous  is  possible.  The  aim  point  for  direct  ascent  during  this  period  of 
time  can  be  continuously  moved  to  that  location  resulting  in  earliest  possible  rendezvous. 

As  the  time  passes  in  such  a  situation  the  rendezvous  point  moves  down- range  and  event¬ 
ually  passes  a  point  where  accomplishment  of  earliest  rendezvous  requires  that  additional 
burns  must  be  added. 

The  approach  is  to  develop  the  program  for  minimum  time  trajectories  for  direct  ascent 
to  rendezvous  for  a  fixed  propellant  budget.  In  addition,  part  of  the  work  of  this  task  will 
be  to  develop  the  policy  for  transition  from  the  direct  ascent  to  the  more  general  variable 
point  policy. 

The  program  would  be  exercised  by  simulation  of  particular  cases  and  the  results  would  be 
evaluated. 

f.  Reentry  Guidance  Implementation  .  A  tie  in  with  re-entry  guidance  should  be  in¬ 
cluded  to  provide  maximum  flexibility  in  mission  application. 

The  descent  problem  may  be  divided  into  two  distinct  steps  that  should  be  considered: 

1. )  Defining  the  pre- entry  trajectory 

2.  )  Guidance  through  the  atmosphere  to  a  landing  point. 

The  primary  consideration  in  planning  the  pre-entry  trajectory  is  the  entry  conditions  in 
velocity  and  flight  path  angie  that  will  be  encour  '  -,red  by  the  vehicle  at  entry  to  the 
atmosphere. 
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Of  special  consideration  in  this  work  is  the  requirement  for  quick  return  of  the  vehicle 
to  a  specified  landing  site. 

The  results  of  the  Variable  Point  Guidance  and  Targeting  additions  necessary  for  pre¬ 
entry  trajectory  determination  will  be  to: 

1.  )  Guide  the  vehicle  to  a  safe  entry  condition 

2.  )  Position  the  nominal  entry  plane  so  that  the  vehicle  can  use  its 

lateral  ranging  capability  to  insure  landing  at  the  specified  site. 

g.  Displays  and  Manual  Control.  Displays  for  both  launch  operations  and  manned 
vehicles  can  be  developed  from  the  basic  Variaole  Point  Guidance  and  Targeting  system. 
The  purpose  of  this  task  is  to  define  the  information  content  and  gross  form  of  the  displays 
required  to  allow  both  the  astronaut  and/or  the  launch  conductor  to  monitor,  modify  or 
control  system  operation  and  performance. 

h.  Error  Analysis.  Guidance  and  Targeting  Error  Analysis.  The  purpose  of  this 
task  is  to  provide  a  method  for  evaluating  the  accuracy  with  which  the  system  is  able  to 
meet  the  desired  mission  objectives.  The  following  subtasks  should  be  performed: 

1)  Navigation  Error  Source  Model.  A  mathematical  model  will  be  developed  to 
represent  the  errors  originating  in  the  navigation  system.  This  will  include  inertial  in¬ 
strument  errors  in  the  IMU,  meaurement  errors  in  the  auxiliary  sensors  used  for 
correcting  the  IMU,  errors  in  ground-a,  rived  data  on  vehicle  and  target  orbit  elements, 
and  errors  due  to  approximations  in  the  navigation  equations,  and  round-off  and  truncation 
in  the  navigation  program. 

2)  Guidance  and  Targeting  Error  Source  Model.  A  mathematical  model  will  be 
developed  to  represent  tne  errors  originating  in  the  guidance  and  targeting  system.  This 
will  include  approximations  in  the  prediction  of  the  future  vehicle  and  target  states,  and 
truncation  and  roundoff  ei  rors  in  the  targeting  and  guidance  programs. 

3)  Control  Error  Sources  Model.  A  mathematical  model  will  be  developed  to 
represent  the  errors  originating  in  the  control  system  for  each  stage.  This  will  include 
uncertainty  of  thrust  tail-off  control  errors  remaining  at  engine  cutoff,  resulting  from 
control  system  lags,  aerodynamic  parameter  uncertainties  and  errors  in  the  atmospheric 
model  as  they  affect  aerodynamic  control. 

4)  Error  Propagation  Model.  A  mathematical  model  will  be  developed  to  represent 
the  propagation  of  the  errors  due  to  the  above  sources  throughout  the  mission.  This  will 
include  both  deterministic  and  stochastic  errors  and  will  present  the  resultant  overall 
error  in  terms  of  deviation  from  the  desired  state  at  critical  mission  times. 
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APPENDIX  I 


DERIVATION  OF  EQUATIONS  OF  ORBITAL  MOTION 


1.  Introduction.  This  appendix  will  derive  certain  equations  pertaining  to  the  motion  of  a 
particle  in  an  inverse  square  centered  force  field.  The  form  of  the  equations  is  chosen  to 
conveniently  fit  the  Variable  Point  Guidance  and  Targeting  policy.  Many  of  these  equations 
arc  standard  in  celestial  mechanics;  a  review  for  the  interested  reader.  Other  equations 
appearing  here  are  not  found  in  the  general  literature  on  orbital  mechanics,  but  have  been 
derived  specifically  to  simplify  the  programming  of  real-time  space  guidance  systems. 

Consider  a  central  force  field  in  which  a  particle  is  at  a  distance  r.  The  acceleration  of 
this  particle  can  then  be  expressed  as: 


R  =  f (R)  Tr 


(AI-1) 


where: 


t(R)  1R  in  any  force  directed  along  R. 


Taking  cross  product  of  both  sides  of  Equation  (AI-1)  with  R  gives 


R  x  R  =  0 


(AI-2) 


which  is  equal  to: 


dF  '  R  x  R  )  =  0 

Integrating  Equation  (AI-3)  gives: 


(Ai-3) 


(AI-4) 


121 


where  h  represents  f.  constant  vector  and  is  the  angular  momentum  of  the  orbit  of  the 
particle. 

Since  h  is  a  constant  vector  and  is  perpendicular  to  the  orbit  ( Ti  and  R)  the  orbit  defines 
a  plane  in  space  which  passes  through  the  force  source. 

The  definition  of  the  cross  product  gives: 

I'RilRjsin  3  =  I  hi  =  Rv  cos  y  =  R2  9  (AI-5) 


where  3  is  the  angle  between  the  radius  and  velocity  vectors.  The  angle  y  represents 
{90  -  3)  and  is  measured  from  the  local  horizontal  to  the  velocity  vector. 

Equation  (AI-5)  can  be  solved  for  time  if  R  is  known  as  a  function  of  0  (equation  of  the 
path) 


dt 


R 


R  V  cos  y 
oo  '  o 


d  0 


(AI-6) 


The  path  equation  can  be  defined  if  the  force  field  characteristics  are  known.  For  the  in 
verse  square  central  force  field  the  acceleration  is 


*  =  '  -T  IT  <AI-7) 

where  jj.  is  the  product  GM  (product  of  gravity  constant  and  mass  of  planet) 

From  Equation  (AI-4)  R  x  R  =  h  which  can  be  crossed  with  both  sides  of  Equation  (AI-7) 
to  give 


R 


x  h 


\  X  (R  x  R) 
RJ 


-Aj-  ( R  .  R)  R  -  (R  .  R)  R 
RJ 


Equation  (AI-8)  is  ther  simplified  giving 


(aI-8) 


R  x  h 

M 


R_  _R  R 
R  '  R  R 


(AI 


since  Equation  (AI-9)  is  known  to  be  the  following  derivative 

/-  -\  /-\ 

d  |R  x  h  |  _  d  /  R  ] 

dt  \  M  /  ~  dt  \R/ 


lAI-10) 


Equation  (AI-9)  may  be  integrated  to  give 


R  x  h 
M 


P.  T 
—  +  e  I 


R 


(AI-11) 


where  e  is  a  constant  and  lc  is  a  unit  vector  along  e.  Multiplying  Equation  (AI-11)  by  R. 


R  .  R  x  h  R  R 


M 


+  e  R  ,  1 
R  e 


(AI-12) 


also 


R.Rxh  RxR.h  h.h 


M 


(AI-13) 


2 

hVn  is  a  new  constant  which  may  be  defined  as 


P 


(AI-14) 


Using  Equations  (AI-13),  (AI-14)  and  simple  vector  algebra,  we  can  rewrite  the  vector 
Equation  (AI-12)  as  the  following  scalar  equation: 

P  -  R  +  eRcos  v  (AI-15) 

where  is  the  angle  between  T  and  5 .  Solving  for  R,  we  obtain 

G 


R  = 


_ P _ 

1  +  e  cos  v 


'AI-16) 


2.  Orbit  Equations.  By  comparing  Equation  (AI-16)  with  the  equation  of  a  general  conic 
section  in  polar  coordinates,  we  identify  e  as  the  eccentricity,  u,  as  the  true  anomaly, 
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and  P  as  the  semi-hitus  rectum  of  the  conic  section  involved.  Substituting  the  expression 
for  "h"  from  Equation  (AI-5)  into  Equation  (AI-14),  we  obtain  another  useful  form  for  P  as 
follows: 


P  = 


d2  2 
R  \T  cos  y 

M 


R 


A 

V  cos  y 

JL 

R 


(AI-17) 


It  will  be  noted  that  ii/R  is  the  square  of  the  circular  velocity  at  a  radius  R.  Denote  the 
circular  velocity  at  radius  Rby  Vo  With  this  definition,  division  of  Equation  (AI-17) 
through  by  R  leads  to 


JP 

R 


2 

cos  y 


?  9 

u  cos  y 


(AI-18) 


where  u  =  V/Vc.  From  Equation  (AI-18)  the  horizontal  component  of  velocity  is  obtained 
from 


(AI-19) 


For  any  point  on  the  conic,  such  as  an  aim  point,  Equation  (AI-16)  may  be  written 


R  = 
a 


1  +  e  cos  v 


_ P 

1  +  e  cos  (v  +  0) 


(AI-20) 


where  v  is  the  true  anomaly  to  the  aim  point  and  Q  is  the  downrange  angle  between  the 
present  point  (r,  v)  and  the  aim  point  (Ra,  ya). 

Expanding  Equation  (AI-19)  and  converting  to  non-dimensional  form 


V_ 

_ R _ 

1  +  e  cos  v  cos  Q  -  e  sin  v  sin  Q 


(AI-21) 


3.  Path  Equations  in  Terms  of  End  Conditions.  In  order  to  write  Equation  (AI-21)  in 
terms  of  initial  conditions  only,  e  cos  v  and  e  sin  v  should  be  written  in  terms  of  initial 
conditions.  The  downrange  angle  ft  may  be  obtained  from  the  dot  product  to  present 
position  and  aim  point  vectors.  From  Equation  (AI-16) 
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e  cos  v 


P 

R 


1 


(AI-22) 


The  function  e  sin  v  is  obtained  from  the  derivative  of  Equation  (AI-16) 


dR 

Rcfy 


tan  y 


e  sin  v 
T~+  e  cos  v 


(AI-23) 


and  since 


P 

R 


1  +  e  cos  v 


e  sin  v  -  tan  y 

Substituting  Equation  (AI-22)  and  (AI-23)  into  Equation  (AI-21)  gives 


(AI-24) 


P 

R 

\  p 

-  1)  cos  ft  -  -p  tany  sin  0 


(AI-25) 


from  Equation  (AI-25)  the  semi-latus  section  ratio,  or  the  horizontal  velocity  ratio  squared 
is  obtained  in  terms  of  initial  conditions. 


P 

R 


1  -  cos  ft _ 

cos  ft  +  sin  ft  tan  y 


and  similarly 


P  1  -  cos  ft 

R  =  R 

-  cos  ft  -  sin  ft  tan  y 

K  a 


(AI-26) 


Equating  P  in  the  above  two  Equations  (AI-26)  gives  the  relation  between  any  two  flight  path 
angles  on  a  trajectory. 
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(AI-27) 


tan  y 


sin 

1  -  cos  Q 


tan  y 
'a 


4.  Required  Velocity  at  the  End  Points.  The  horizontal  component  of  velocity  at  the  point 
at  radius  R  was  given  in  Equation 19)  as 


V„  =  V  cos  y 


1/2 


The  radial  component  is  then  obtained  by  combining  Equation  (AI-19)  and  (AI-27)  which 
gives  the  flight  path  angle  relationship  at  the  end  points 


VR  =  V  sin  y  =  tan  y 


(AI-28) 


The  velocity  components  entering  the  final  point  are  then: 


VTT  =  V  cos  y  = 
Ha  a  'a 


_Z_  JL 

R  R  j 
\  a  a  / 


y/2 


R 

R 


V  cos  y 


and 


VRa  =  Va  sin  va  =  VHa  tan  ra 


(AI-29) 


The  required  reference  velocity  components  which  are  used  to  derive  Vg  are 


l/2 

V  cos  y  =(-£•#)  =  V„ 

req  req  l  R  R  j  H  req 


and 


V  sin  y  =  VXT  tan  y 
req  req  H  req  req 


(AI-30) 


transforming  to  earth  centered  coordinates 


V  =  1D  V  sin  y  +  1„  V  cos  y 
req  R  req  req  H  req  'req 


(AI-31) 
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u-atm t#******- - 


where  Tr  and  are  unit  vectors  along  the  radius  vector  and  perpendicular  to  the  radius 
vector  in  the  guidance  plane  or 


XH  *h  X  *R 


>h  = 


TRxIRa 

iTrxTr1 


(AI-2?) 


Vg  is  then  derived  by  subtracting  the  components  of  present  velocity  obtained  from  the 
navigation  computation. 


V  =  V  -  V  (AI-33) 

g  req 


3.  Tangent  Transfer  Trajectories.  Equation  (AI-27)  above  can  be  rearranged  to  give: 


tan  y  -  -w-  tan  y 

sin  ft  _ _ _a _ 

1  -  cos  ft  ~  -  R 

’  R 

a 


(AI-34) 


or 


tan 


ft 

2 


tan  y  -  tan  y 
Ra 


a 


Equations  (AI-34)  are  useful  for  deriving  the  tangent  transfer  trajectory  equations.  The 
end  points  and  the  included  angle  ft  for  the  best  trajectory  need  to  be  found  without 
ambiguity  or  quadrant  problems.  Analysis  and  experiment  have  shown  that  the  tangent 
transfer  trajectory  between  two  orbits  is  efficient  in  propellant  usage  if  the  tangency  at  the 
higher  altitude  end  is  placed  at  a  node.  Also  the  most  efficient  node  is,  in  general,  the 
higher  of  the  two  possible  nodal  points.  Results  of  detailed  studies  on  the  efficiency  of 
the  tangent  transfer  were  given  in  reference  _9 _ . 

For  circular  orbits  or  for  elliptical  orbits  whose  major  axes  are  both  aligned  with  the  line 
of  nodes  the  tangent  transfer  from  node  to  node  with  optimal  plane  change  split  gives  the 
optimal  solution.  If  the  major  axes  are  not  aligned  or  if  one  ellipse  is  not  aligned  with  the 
line  of  nodes,  the  tangent  transfer  at  the  highest  node  is  c*ose  to  the  optimum  transfer 
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and  the  determination  of  the  end  points  for  the  general  case  is  given  as  follows: 

Figure  58  illustrates  the  general  case  of  the  tangent  transfer  situation  when  both  orbits 
are  elliptical  and  not  aligned. 


Figure  61.  Tangent  Transfer  Trajectories 
Equation  (AI-27)  which  gives  the  downrange  angle  function 


tan 


JL 

Z 


is  rewritten  in  terms  of  the  nomenclature  given  in  Figure  61. 


1^ 

Z 


RNt 

R. 

1 


tan  y 


Nt 


-  1 
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l 


tany. 


Having  defined  the  tangent  transfer,  it  follows  that  any  tangent  trajectory  can  readily  be 
defined  traversing  from  or  to  any  initial  or  final  point  on  either  trajectory  (the  transfer 
always  exists  for  non -intersecting  orbits.  For  intersecting  O  'bits  one  of  the  nodes  will 
produce  a  workable  solution).  Since  the  angle  from  any  point  on  one  of  the  orbits  be¬ 
tween  the  lines  to  the  empty  foci  is  equal  to  twice  the  flight  path  angle  at  the  point,  the 
other  end  point  for  the  tangent  ti  ijectory  is  easily  constructed. 

The  outer  orbit  in  Figure  58  is  considered  to  be  the  target  trajectory  with  no  particular 
time  or  position  specified.  F’  is  the  location  of  its  empty  focus.  F  is  the  center  of  the 
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force  field  and  is  thus  one  focus  for  any  orbit  in  this  field.  The  inner  orbit  is  the  inter¬ 
ceptor  trajectory  with  empty  focus  at  Fj.  and  are  radius  vectors  along  the  line  of 
nodes  between  the  two  orbits.  R^ji  is  the  nodal  interceptor  radius  vector  and  R^t  Is  the 
nodal  target  radius  vector. 

Rt  and  Rj  are  target  and  interceptor  radius  vectors  at  an  angle  ft  from  the  line  of  nodes. 
Rt  is  in  the  target  orbit  plane  and  Rj  is  in  the  interceptor  orbit  plane. 

Four  trajectories  can  be  passed  between  the  four  points  shown  above.  The  object  of  the 
proof  here  is  to  show  that  all  four  of  these  trajectories  contain  the  same  angle  ft  and  are 
all  tangent  to  each  other  at  the  indicated  end  points.  We  have  already  defined  one  of  these 
trajectories. 

Two  other  tangent  trajectories  are  tangent  by  definition.  They  both  lie  on  the  two  original 
orbits  and  in  their  own  respective  planes.  The  same  angle  is  traversed,  from  a  crossing 
vector  (the  line  of  nodes)  to  another  point  along  their  original  orbits.  The  above  three 
trajectories  contain  the  same  central  angle  and  have  end  points  of  radi  as  stated. 

As  previously  stated  the  fourth  trajectory  is  considered  as  one  starting  in  the  target  plane 
at  point  Rfc  and  arriving  at  R^i  in  the  target  orbit  plane.  To  prove  that  these  four  tra¬ 
jectories  have  the  same  flight  path  angles  at  their  respective  end  points  which  contain  the 
same  included  angle  requires  that  the  tangency  conditions  exist  for  the  fourth  specific  tra¬ 
jectory  just  defined. 

Let  X  =  tan  y  and  write  Equation  (AI-28)  for  each  trajectory 


R. 

1  Y 
RXT.  Ni 
Ni 


also 


X. 

1 


1  •  ir- 

RNt 


x 

RNt  Nt 


(AI-35) 


and 
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For  the  fourth  trajectory 


R. 


1  - 


t 


R 


_ Ni 

(tan 


R 


R 


Ni 


(A.I-C6) 


To  prove  the  conjecture  equate  the  first  of  the  above  two  equations  and  solve  for  then 
substitute  the  solution  fcr  Xjvjt  from  the  third  into  the  last  equation.  This  gives  the  desired 
result 


Further  algebraic  manipulation  of  these  equations  then  gives  the  range  angle  ft  as  a  function 
of  the  radii  and  flight  path  angles  of  each  orbit  along  a  node  as  follows: 


Tan 


ft 


«Ni 


-  1 


^Nt 

tan  rNt  -  R7T.  tan  yNi 
Ni 


(AI-37) 


since: 


.  0 
tan  T 


1 

z 


sin  ft 


2  Z 


12  "  z2  ,  1 


(AI-38) 


cos  ft12  -  1 


Z  +  1 


These  are  convenient  forms  for  use  in  the  program.  The  details  of  application  are  given 
in  Volume  I  part  2. 

6.  Time  of  Flight.  Up  to  this  point  the  equation  of  a  path  has  been  specified  with 
assoc iatecf initial  conditions  which  pass  through  a  specific  point  in  space,  also  the  tangent 
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C/ 


transfer  trajectory  parameters  have  been  derived  but  one  or  more  paths  may  be  required 
which  fits  an  overall  time  constraint. 

In  solving  Equation  (AI-6)  for  time,  using  the  path  equation,  it  is  useful  to  first  derive 
some  of  the  orbit  elements  in  terms  of  end  conditions. 

From  Equation  (AI-22)  and  (AI-24) 


2  2  2 
e  =  (e  sin  v)  +  (e  cos  v) 


or 


2  P 
e  =  It 


(1  +  tan2  y)  -  2 


+  1 


also  from  Equation  (AI-18)  and  (AI-28) 


2  -  2  2  (r  2X 
e  =  1  -  u  cos  y  (2  -  u  ) 


(AI-38) 


(AI-39) 


The  major  axis  ratio  is  then  found  from  Equation  (AI-16)  as  the  mean  of  perigee  and 
apogee  distance 


a_ 

R 


(AMO) 


in  terms  of  initial  velocity 


_ 1 

2  -  ~  (1  +  tan^  y) 


(AI-41) 


from  Equation  (AI-39) 


2 

e  = 


(AI-42) 


131 


and 


b2  =  a2  (1  -  e2)  (AI-43) 

then 

(|)2  =  (?)  (AI-44) 

also 

( <«-«) 

The  time  of  flight  from  start  point  to  aim  point  can  now  be  found  by  integration  of 
Equation  (AI-6). 

First  let 


T  = 


R 


(AI-46) 


for  non-dimensional  time.  Then  equation  (AI-6)  in  non -dimensional  form  is: 


(AI-47) 


Then  from  Equation  (AI-20)  and  (AI-25) 


T 


lA2/2t 

d  0 

\rJ  Jo 

1  ♦(£-*) 

_  j 

cos  0  -  -5  tan  y  sin  0 

rt 

(AI-4P) 


/p\3/V°o+fi 
\R)  •(>  1  + 


d  (o  +  eo) 


e  cos  (O+O) 
o 


(AI-49) 
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from  Pierce’s  tables  with 


e  <  1 


(AI-50) 


The  time  equation  then  may  be  written  for  any  start  point  to  any  angle  to  2 n  as 


*  -(*) 


3/2  r 


2  arc  cot  i  (l  !ircofls  o  - tan  y) 


-  S  (tany  -  tany) 

a  a 


(A  1-51) 


Similar  solutions  of  the  time  equation  are  readily  obtained  for  the  cases  where  e  =  1  and 
e  >  1.  For  the  case  where  e  >  1 


>  i  =  (i) '  [- 2  arc  coth  1  (r-~“ 03  a-  - tan  y) 


+  |  (tan  ya  -  tan  y) 


(AI-52) 


which  gives  the  time  of  flight  on  an  hyperbolic  trajectory  and  for  e  =  1 


e  =  1  3  \  R 


i/p\1/2[p.a  /  p 

-  (tan  yal  ^1  +  p- 


W 

[“R 


-  (tany)  U  +  r- 


(AI-53) 


which  gives  the  time  of  flight  on  a  parabolic  trajectory.  This  equation  is  useful  when  the 
ratio  R/a  approaches  zero. 

Since  for  the  parabolic  trajectory  the  velocity  anywhere  along  its  path  is  equal  to  the  escape 
velocity,  the  seir.i-latus  rectum  ratio  is  a  function  of  the  radius  vector  to  the  point  and  the 
flight  path  angle  at  the  point  or 


JP 

R 


=  2 
e  =  1 


2 

cos  y 


: — — r- 

1  +  tan  y 
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and 


„  2 

2  cos  y  -  • 

T  2L 

1 

In  terms  of  flight  path  angles, 


2 

.  ,  2 
1  +  tan  y 

a 

range  angle  and  radius  ratio 
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tan  y 


where 


tan  y  = 
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sin  Q  /Ra 
1  -  cos  ft  \  R 


R 

if  tany 


(AI-54) 


and 


7.  Time-To-Start-Burn  for  a  Two -Stage  Vehicle.  One  of  the  basic  functions  of  Variable 
Point  Guidance  is  to  determine  the  remainder  of  an  ascent  trajectory  based  on  present 
situation  parameters.  This  trajectory,  if  successfully  executed  by  the  vehicle  in  response 
to  the  commands  of  the  guidance  computer,  will  lead  to  rendezvous  with  the  intended 
orbital  target.  Determination  of  an  ascent  trajectory  and  issuance  of  suitable  commands 
to  achieve  this  trajectory  are  continually  performed  by  the  Variable  Point  Guidance  Com¬ 
puter  while  the  vehicle  is  both  on  the  pad  and  in  its  ascent  flight.  The  ascent  trajectory 
decided  upon  is  computed  based  on  impulsive  burns  located  at  aim  points  in  space  which 
are  varied  in  position  to  satisfy  the  time  constraint  for  rendezvous  (within  certain  vehicle 
safety  constraints).  The  number  of  aim  points  required  is  determined  by  the  guidance 
computer. 

As  the  intercepting  vehicle  approaches  an  aim  point,  the  guidance  computer  must  decide 
when  to  start  the  (finite)  vehicle  thrust  to  match  the  effect  of  the  impulsive  maneuver 
planned  to  take  place  at  that  aim  point.  It  can  be  shown  that  under  generally  applicable 
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conditions,  a  finite  burn  can  be  made  equivalent  to  an  impulsive  maneuver  (requiring  zero 
time)  if  the  finite  burn  is  begun  at  the  proper  time  and  is  properly  oriented  through  its 
duration. 

In  this  section  the  correct  time  to  start  such  a  burn  is  derived  for  the  case  in  which  the 
maneuver  is  to  be  executed  by  a  two- stage  rocket  vehicle,  each  stage  of  which  has  a  pre¬ 
dictable  initial  thrust,  a  thrust  time-history  predictable  by  the  rocket  equation,  and  known 
cutoff  conditions. 

Consider  the  vehicle  as  initially  on  a  path  which  if  continued  unaltered  would  pass  through 
the  computed  aim  point  at  the  correct  time  associated  with  that  aim  point.  The  trajectory 
leaving  this  point  must  then  be  changed  to  one  that  will  arrive  at  the  next  selected  aim 
point  at  its  assigned  time.  Figure  59  illustrates  the  situation  at  this  time.  All  vectors 
shown  represent  3x1  column  matrices  in  an  earth-centered  inertial  coordinate  system. 
As  the  vehicle  approaches  aim  point  "a”  the  guidance  computer  attacks  the  problem  of 
reaching  aim  point  "a  +  1"  on  time.  The  known  conditions  are: 

1. )  Rq  =  present  vehicle  position 

2. )  Vo  =  present  vehicle  velocity 

3. )  Tq  =  present  time 

4. )  Ra+  j  =  aim  point  location 

5. )  Ta+  j  =  time  at  aim  point  "a  +  1” 

6. )  ya+1  =  flight  path  angle  at  arrival  at  "a  +  1”. 

From  the  above,  the  guidance  computer  determines  Vreq,  the  velocity  necessary  at  the 
present  time  to  arrive  at  ”a  +  1”  with  a  flight  path  angle  of  ya+l-  The  velocity  vector  that 
must  be  added,  Vg,  is  therefore  computed  by: 


(AI-55) 


If  the  vehicle  thrust  acceleration  amplitude  profile  is  a  known  function  of  time,  a(t),  and  the 
thrust  acceleration  vector  is  kept  parallel  to  Vg,  the  predicted  time  required  for  thrusting, 
T,  is  defined  by  the  scalar  equation: 


(AI-56) 


The  position  of  the  vehicle,  Rco,  at  the  predicted  end  of  the  thrust  period  is  given  by: 
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Figure  62.  Time  to  Start  Burn  Situation  Parameters 
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(AI-57) 


R  =  R  +  V  T  + 
co  o  o 


foil  i(t)dt  Em+  E?(t) 


where: 


:r<‘>  \f 0  /  0 


GM  ,  / V(t)  •  R(t)\2  1 


0  l  R(t) 


J  +  V  R(tr  /  Rft) 


dt  dt 


and 


E77 


rr  r{  L  v(t)  •  R(t)  \ 


v(t)=/0/0\vo-Ttf^; 


dt  dt 


The  effects  introduced  in  the  cutoff  position  (R^)  by  errors  E^(t)  along  R(t)  and  Eym 
along  V(t),  which  are  due  to  gravity  during  the  burn,  are  not  considered.  This  approxi¬ 
mation  is  justified  to  the  extent  that  the  burn  duration  is  short  and  the  gravity  gradient 
associated  with  the  difference  in  positions  at  initiation  and  cutoff  is  small. 

Equation  (AI-57)  may  then  be  rewritten  as: 


a(t)  dt  dt 


(AI-58) 


If  the  time  interval  from  start  of  thrust  to  the  time  the  impulse  was  to  occur  is  called  T  , 
the  desired  position,  R^,  at  the  end  of  time  interval  T  (i.e. ,  that  where  an  impulsive  S 
vehicle  would  be)  is: 


Rd  (T)  =  Rq  +  VQ  Tg  +  VT  (T  -  Tg)  (AI-59) 


The  time  to  start  the  burn  is  then  found  by  equating  Equation  (AI-58)  and  (AI-59),  setting 

R  equal  to  R  .  (T). 
co  d 


Rd(T>  =  Rco 

R„  +  V  T  +  V_,  (T  -  T  )  =  R  +  V  T  +  [  [  a(t)  dt  dt  fAI-60) 

cost  s  0  0  J  Oj  0 
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Equation  (AI-60)  is  illustrated  in  Figure  60  in  which  both  the  impulsive  vehicle  and  the 
real  vehicle  are  shown  in  the  same  positions  with  the  same  velocity  at  the  times  of 
initiation  and  cutoff  of  the  burn.  The  components  of  distance  which  comprise  the  total 
position  change  are  also  shown.  Equation  (AI-60)  reduces  to: 


t=  =  t  ■  t/o71 


a(t)  dt  dt 


(AI-61) 


Interval  Ts  is  the  ’’lead  time”  required  in  the  start  of  the  burn,  T  being  computed  from 
Equation  (AI-56). 

A  two-stage  vehicle  with  constant  thrust  and  mass  flow  rate  for  each  stage  is  assumed.  The 
thrust  parameters  for  each  stage  are  represented  as  follows: 

Stage  1  Stage  2 

Initial  Thrust  Acceleration  a  -  a  n 


Exhaust  Velocity  Cj  C ^ 

Burn  Time  Remaining  Tj  Tg 

If  the  second  tage  is  started  at  some  time  aT  after  the  remaining  burn  time,  Tj,  of  the 
first  stage  the  vehicle  acceleration  profile  at  a  time  t  (measured  from  initiation)  will  be: 


t  <  T« 


(AI-62) 


a  =  0, 


Tj  t  Tj  +  aT 


(AI-63) 


,  Tj  +  ^T  <  t  <  T 


(AI-64) 


ft’ VAT) 


By  integrating  the  acceleration  profile,  the  velocity  gained  during  use  of  the  first  stage 
becomes: 


v,  =  -  c.  In  [l  -  ^7]  ,  t  <  T. 


(AI-65) 
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Figure  63*  Equivalent  Position  Burn  Computation 
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No  velocity  is  gained  through  the  "dead  time",  aT  (referencing  the  development  of  Equation 
(Al-58).  The  velocity  gained  through  use  of  both  the  first  and  second  stages  is: 


VT  =  V1  +  V2 


VT  =  V.,  +  C2  In 


1  -  cf 


,  T-  +  AT  <  t  <T 


(AI-66) 


Equation  (AI-62)  n.ay  be  used  to  find  thrust  time  T  in  terms  of  Vt  as: 


T  = 


ol 


— ® 


,  v  <  V- 
’  g  i 


(AI-67) 


T  =  T.  +  AT  +  - 
1  a 


o2 


1  -  exp  (  - 

u2 


V  >  V, 


(AI-68) 


Integrating  Equation  (AI-65)  to  determine  the  distance  S  covered  during  the  burns  as  a  re¬ 
sult  of  burn  acceleration  only  gives: 


S  =  f  f  a(t)  dt  dt  =  f  V(t)  dt 

J  0  J  0  Jo 


(AI-69) 


S  = 


*ol 


a  i  t 

i  j  ol 
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a  -  t 

In  (1  -  —  1  -  1 
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t  <  T 

’  1 


(AI-70) 


S  =  S, 


S  =  Sj  +  Vj  (t  -  Tj) 


t  =  T, 


Tj  <t  <  Tj  +  AT 


(AI-71) 

(AI-72) 
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(AI-73) 


where: 


tj  =  At  <t  <  T 


Expressing  S  in  terms  of  T  and  substituting  into  Equation  (AI-61)  yields: 


(AI-74) 


(AI-75) 


It  can  be  seen  that  the  above  equations  are  readily  adaptable  to  computer  solution.  A 
possible  computation  sequence  for  the  situation  shown  in  Figure  1  is  as  follows: 

fr>  From  *o>  V  To>  Ta+1’  \+V  awlJ'a+l’  compile  VREQ  and  Vg 

(b)  Compute  present  stage  velocity  capability  remaining,  V.,  using 
Equation  (AI-65) 

(c)  Compare  to 

1. )  If  V1  >  Vr,  compute  time  to  start  burn  fT  )  using  Equation  (AI-74). 

i  g  s 

2. )  If  V.,  <  V  ,  compute  T  using  Equations  (AMO)  and  (AI-75). 

i  g  s 

(d)  Compare  Tq  to  (Ta  -  Tg) 
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1.) 

2.) 


If  T  <  (T  -  T  ),  do  not  thrust 
o  a  s 

H  To  >  (Ta  -  Tg),  start  engine; 
Vg  passes  through  zero. 


until  next  computation  cycle 
cutoff  when  forward  component  of 


After  the  engine  is  started,  the  variable  point  V«  steering  control  will  command  a  \  ''cision 
cutoff  consistent  with  successful  arrival  at  aim  point  "a  +  1". 


APPENDIX  : 

LIST  OF  ABBREVIATIONS  AND  SYMBOLS 


Symbol 

Definition 

A,  B 

= 

Coefficients  depending  on  subscripts 

a 

= 

Semi  major  axis  of  an  orbit 

aft) 

= 

Acceleration  as  a  function  of  time 

b 

= 

Semi  minor  axis  of  an  orbit 

c 

= 

Exhaust  velocity 

E 

= 

Eccentric  anomaly  of  an  orbit 

e 

= 

Eccentricity  of  an  orbit 

F 

= 

Force  vector 

f 

= 

Specific  force  (thrust  to  mass  ratio) 

g 

= 

Gravity  vector 

G 

= 

Gravitational  constant 

H,  h 

= 

Angular  momentum  of  an  orbit 

i 

= 

inclination  of  an  orbit 

W 

= 

Specific  impulse 

J 

= 

Oblateness  coefficient 

K 

= 

Coefficient  or  constant  depending  upon  subscript 

L 

Lattitudc  of  planet  surface  location 

M 

= 

Mean  anomaly,  mass  of  planet 

m 

Mass  of  vehicle 

n 

= 

Number  of  orbits  (parking  or  phasing) 
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5 


Symbol 

P 

R 

T 

u 

V 

V 

req 


Definition 

Semi  latus  rectum  of  an  orbit 

Radius  vector  from  center  of  force  field 

Time  or  time  interval  denoted  by  subscript 

Velocity  ratio 

Velocity  vector 

Velocity  required 

Velocity  to  be  gained 


a 

0 

y 

6,  A 


A 

P 

v 

P 

T 

U> 

ft 


=  Plane  change  angle 
=  Coefficient 
=  Flight  path  angle 
=  Finite  difference 
-  Lagrange  multiplier 

=  GM  =  gravitational  constant  for  specific  planet 
=  True  anomaly  of  point  on  an  orbit 
=  Relative  range 
=  Period  of  an  orbit 
=  Argument  of  perigee  of  an  orbit 

=  Angle  between  two  points  on  an  orbit,  also  the  right  ascention  of  the  ascending 
mode 


Subscripts 
at  = 

a  = 

B 
H 
ij 


Refers  to  apogee 
Refers  to  aim  point 
Body  axes 
Horizontal 

Points  at  beginning  and  end  of  arc 
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Symbol 

P 

M 

7 T 

r 

t 

<t> 


x,y,z 


4,  v>  £ 


Definition 

Refers  to  platform  axes 
Refers  to  intercepting  vehicle 
Refers  to  perigee 
Point  R 

Refers  to  target  vehicle 
Refers  to  phasing  unit 
Earth  centered  inertial  (ECI)  axes 
Steering  axes 


COORDINATE  SYSTEMS 


XYZ  ECI  COORDINATES 
XB  YB  ZB  BODY  COORDINATES 

Figure  64.  Coordinates 


Figure  65,  Orbit  Parameters 
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